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ABSTRACT 


In order to study the nucleate pool-boiling performance 
of R-114 and R-114-o0il mixtures from enhanced evaporator 
tube surfaces, an experimental apparatus was designed, 
constructed and instrumented. The evaporator was made of a 
T-shaped Pyrex glass container. Boiling occurred from a 
smooth, hard-copper tube, 15.9 mm (5/8 in.) in outer diam- 


eter, 12.7 mm (1/2 in.) in inside diameter and 43126 mnie 


Teo in length. The tube was heated uSing a cartridge 
heater, and it was instrumented with 8 thermocouples to 
measure the wall temperature. A Hewlett-Packard 3497A data 


acquisition/control unit and a 9s826A CGomputer were  uUsed@ee 
collect and process data. The condenser was cooled by an 
ethylene glycol-water mixture, which was maintained at about 
-17 °C by means of an R-12 refrigeration system. Nine data 
runs were completed to de-bug the experimental apparatus and 
to check for réeprodtier bite. During all data runs, espe. 
cially at higher heat fluxes (greater than 10 kW/m?), large 
temperature variations were observed along and around the 
active boiling length of the test tube. The data were 
compared with data found in the literature and reasonable 


agreement was obtained. 
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le 2NTRODUGT ON 


A. BACKGROUND 


In order to reduce the size and weight of “ailing 
conditioning systems and water chillers, the United States 
Navy recently has considered two different approaches: 1) to 
change the type of the refrigerant from R-1ll to R-114, and 
2) to use enhanced evaporator tube surfaces, such as Gewa-T, 
Thermoexcel-E, and High-Flux. Schematic drawings and photo- 
micrographs of these enhanced surfaces at various magnifica- 
tions are given in Figures 1.1 ,1.2 and 1.3. <A combingedenw 
of R-114 and enhanced surfaces may lead to significant 
improvements in system efficiency as well as to smaller and 


lighter components in future Navy applications. 


B. ADVANTAGES OF USING R-114 


The latent heat per unit volume of R-114 is 50 BTU/ft? 
at 70 °F. As a result, the heat carried by each cubic foot 
of R-114 is nearly twice as much as that carried by R-l1l 
[Ref. l] at the same temperature. In other wouder the 
volume of R-114 vapor circulated per ton of refrigeration is 
about half that for R-ll. Inherently, the size of the system 
components can be reduced significantly using R-114 rather 
than R-ll. R-12 and R-22 in fact can transfer the heat nome 
efficiently than R-114, but they require higher pressure 
ranges. However, operating pressures for R-114 are moderate, 
Similiar tomas such that light-weight components may be 
desirable and less-experienced personnel may be involved in 
the installation and maintenance. While the refrigerants in 
general are nonflammable, R-114 and R-12 are rated to have 


the lowest toxicity among all refrigerants [Ref. 2]. Based 
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on the above discussion, R-114 is probably the most suitable 
refrigerant for use on board submarines and ships. 

Further, R-114 is considerably more stable toward heat 
mean R-1l2. All standard ferrous and nonferrous’ metals, 
except zinc and magnesium, and aluminum alloys containing 
appreciable amounts of zine and magnesium, may be used with 
PemelG |Ref. 3]. The physical properties of R-114 and its 
comparison with the other common refrigerants are summarized 
in Appendix A, pressure-enthalpy diagram of R-114 is 
provided in Appendix B andthe range of application of 


refrigerants are provided in Appendix C [Ref. 4]. 


C. ENHANCED BOILING TUBE SURFACES 


The evalution of enhanced Surface geometries that 


promote high-performance nucleate boiling has been surveyed 


by Webb [Ref. 5]. Two basic types of enhanced boiling 
surfaces are employed on the outer surface of tubes: porous 
coatings and reentrant grooves. The High-Flux surface is a 


wemeered, porous metallic matrix bonded to the base surface. 
The Gewa-T and Thermoexcel-E surfaces have reentrant-type 
nucleation sites formed by cold working an integral finned 
foe wall. See Figures 1.1, 1.2 and 1.3 for details. 

The key to the high performance of these surfaces is 
attributed to three factors [Ref. 5] : "(1) a pore or reent- 
GEemit cavity within a critical size range, (2) interconnected 
cavities, and (3) nucleation sites of a reentrant shape. 
When the cavities are interconnected, one active cavity can 
activate adjacent cavities. It appears that the dominant 
Meret lon Of the vaporization occurs at avery thin liquid 
film within the subsurface structure. The reentrant cavity 
Shape provides a stable vapor trap, which will remain active 


at very low liquid superheat values.” 


iS 


The pool boiling heat-transfer performance of R-l13  faeg 
the above-mentioned commercially-available enhanced surfaces 
has been investigated by Marto and Lepere [Ref. 6]. They 
reported that for R-l1l3, the enhanced surfaces showed a two 
to tenfold increase in the heat-transfer coefficient when 
compared to a plain tube. The High-Flux surface was most 
effective over a broad range of heat fluxes, whereas the 
Gewa-T surface was only effective at high heat fluxes (near 
100 kW/m?). The Thermoexcel-E surface showed similar gains 
in the heat-transfer coefficient to that of the High-Flux 
surface below 10 kW/m?. They also reported that the degree 
of superheat required to activate the enhanced surfaces was 
less than the plain tube, and was sensitive to initial 


surface jeonditioning Ande tf lutawonope.emece 


D. NUCLEATE POOL BOILING PERFORMANCE OF PURE R-114 


The heat-transfer performance of pure R-114 in the 
evaporator section of an air-conditioning unit Can bewecem. 
mated using one of the Suggested commonly-used correlations 
fous nucleate peel Diode erie which are listed by 
Chongrungreong and Sauer [Ref. 7]. However, because of the 
lack of complete understanding of the nucleate boiling 
mechanism, these correlations should be compared with exper- 
imental data a1f 1t is “poss pales Son i daar the literature 
reveals a limited amount of nucleate boiling experimental 
investigations using R-114. Happel [Ref. 8], using a hori- 
zontal, pure-nickel tube with a surface roughness of 0.43 
micrometer, investigated the boiling characteristics of pure 
R-114 at various pressures. The effect of pressure on 
nucleate boiling heat transfer in R-114 was also investi- 
gated experimentally by Nishikawa and his co-workers 
[Ref. 9]. They operated over “a@uedulécdiip1 es cuiem «i 
range from Q.03 to 0.98, and demonstrated that the heat 
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transfer coefficient of R-114 increases remarkably at high 
pressure near the critical point. Their results showed that 
the effect of surface roughness changes with pressure and 
the effect gradually vanishes as pressure approaches” the 
aioe cal point. Using a horizontal smooth copper’ tube, 
Henrici [Ref. 10] investigated boiling heat-transfer 
performance of R-114 over a heat-flux range from 0.1 to 100 
ilmvea- and over a pressure range from 37 to 252 kN/m?. 
Stephan and Mitrovic [Ref. 11] investigated the perform- 
ance of pure R-114 from a Gewa-T enhanced surface. Boiling 
occurred outside a horizontal bundle of these tubes. They 
reported that the position of the tube within the bundle 
does not effect the heat-transfer coefficient considerably. 
The reason for this effect was explained in the following 
manner: ‘The flow pattern around a tube seems to be of very 
weak influence and the heat-transfer process 1S mainly 
governed by the bubble formation and the two-phase flow 
inside the T-shaped channel. The top part of the fins, which 


were rolled to form. the T, Piefeeitedtivw~encitca =  MOt direct Ly 
contribute to the heat-transfer process becauSe, even at 
maximum heat fluxes, vapour bubbles were not produced 
there." 


Marto and Hernandez [Ref. 12] also investigated nucleate 
meme DOlling Characteristics of a Gewa-T surface in R-113. 
Their experimental results showed a disagreement with 
Seeplian and Mitrovic’s above explanation. They showed that 
“the liquid-vapor motion may not always be confined to flow 
circumferentially within the channels formed by the T-shaped 
fins as previously postulated. Surface imperfections on the 
tips of the Gewa-T fins may cause bubble nucleation from 
these sites. The presence of these active Sites may further 
enhance the performance of the Gewa-T surface." 

The factors which appear to affect boiling heat-transfer 


performance of a Single tube are the physical properties of 
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refrigerant, boiling pressure, tube diameter, surface condi- 
tion of the tube, and maybe the "Submergence effect” 
(resulting from the hydrostatic liquid head) on the evapora- 
tion process. Experimental observations have shown that 
changes in magnitude of these properties and conditions 
Significantly affect pool-boiling heat-transfer performance. 
[Ref. 7]. 


E. NUCLEATE BOILING PERFORMANCE OF REFRIGERANT-OLIL MIXTURES 


Most of the refrigeration systems work WLEh iowa 
JubT Heaedr hermetically sealed compressors, which are 
cooled by the working fluid. Therefore, the fluids being 
circulated within these refrigeration systems generally are 
not pure refrigerants, but are refrigerant-oil mixtures. 
The saturation temperature of a refrigerant-oil mixture is 
higher than that for the pure refrigerant at the same pres- 
sure. Also, both viscosity and surface tension increase 
when oil is added to refrigerants. Other thermodynamic prop- 
erties of refrigerants, such as density and specific heat, 
are also significantly changed by the addition of oi }3) ihm 
adequate knowledge on the influence of oil on the boiling of 
refrigerants 1S very important for the successful design of 
refrigeration systems. The literature reveals only very 
little information on the boiling performance of R-114-o0i1 
MLxbures | Veteemee ln In particular, heat-transfer perform- 
ance of the enhanced evaporator surfaces in refrigerant-oil 
mixtures 1S a poorly-understood Subject area. Since these 
advanced boiling surfaces contain a high density of pores, 
the presence of oil can generally result in a greater 
performance degradation than on more-conventional, smooth 
surfaces. 

Various studies which are surveyed and listed by Jensen 


and Jackman [Ref. 13] have shown that, in general, nucleate 
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boiling heat-transfer coefficients decrease with increasing 
MMMCOncentiaelmonm Lik- ll. R-12, R-22, R-113 and R-114. For 
oil concentrations less than 34, however, Dougherty and 
Sauer [Ref. 14] demonstrated that the heat-tansfer coeffi- 
Beene may be increased slightly for R-li-oil and R-113-011 
mixtures. They attributed this trend to an increase of 
foaming of the mixtures. This result also has been- confirmed 
by Saver, Gibson and Chongrungreong [Ref. 15] using R-12-o0il 


mixtures, and by Henrici [Ref. 10] using R-114-o0il mixtures. 


However, in a recent experimental study using R-1l3-oil 
feeeeeres, Jensen and Jackman [Ref. 13] did not observe such 
an enhancement with low oil concentrations. The model of 


the governing processes involved in boiling of refrigerant- 
Smiemixtures has been given by Stephan [Ref. 16] as follows: 
“Obviously the reduction in heat transfer in a mixture is 
due to the fact that the liquid near the vapour-liquid 
interface of growing bubbles becomes denuded in the more 
volatile components. Thus, a concentration gradient exists 
between the interface and the bulk of the liquid, and the 
more volatile components have to diffuse through the liquid 
in order to reach the interface. A mass-transfer resistance 
thus 1s introduced, which is not present when boiling single 
components.” 

The effect of the oil viscosity on heat transfer in 
boiling of refrigerant-oil mixtures has been considered by 
Samer, Gibson, and Chongrungreong [Ref. 15]. Onli seosa ty. 
1s believed to significantly contribute to the changes in 
Pemeine characteristics of refrigerant-oil mixtures. The 
heat-transfer coefficient 1S greatest Hom eamemm i lLehest 
viscosity of oil for the same temperature difference. 

The effect of oil on boiling of R-12 from a Gewa-T tube 
Was investigated by Stephan and Mitrovic [Ref. 17]. They 
reported that the presence of oil in R-12 reduced the heat- 


transfer coefficient almost over the entire range of heat 


eo 


fluxes and mass fraeGevous. They observed an increase in 
heat-transfer coefficient with less than 54 mass fraction of 
oil only at high heat fluxes (such as higher than 20 vk) ieee 

The effects of oil on the heat-transfer performance of 
refrigerants ina tube bundle have been summarized by Arai, 
eta Le (Ref. 18] in the following manner: "The presence of 
oil in a shell-and-tube evaporator is believed to have two 
opposing effects on the heat-transfer coefficient of the 
evaporator tube bundles. Firstly, many people suspect that 


the oil left out from the evaporating refrigerant tends to 


choke the surface structure, thereby lowering the heat- 
transfer coefficient. Secondly, from the experience of 
conventional heat exchangers, it has long been recognized 


that the presence of oil enhances foaming on the upper part 
of the tube bundles. This foaming increases the vigor of 
refrigerant movement around the tubes, so that heat transfer 


1S ennmemecca:: 


Fo THEsts Chace rr e. 


Based on the foregoing discussion, it 1s clear that 


little information exists in the literature to enab ween 


Successful design of R-114 evaporators, especially when 
uSing enhanced boiling surfaces, in the presence of lubri- 
cating oils. The objectives of this research effort were 


therefore to: 

1. Design, construct and instrument an experimental 
apparatus in order to investigate the heat-transfer 
performance of R-114 and R-114-01i1 mixtures from 
enhanced evaporator tube surfaces, 

2. Successfully operate the system to produce repeatable 
data, and 

3. Obtain smooth-tube baseline data for R-114 to be used 
as a standard for the future data on advanced 


SUriaces. 
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II. DESCRIPTION OF EXPERIMENTAL APPARATUS 


A. OVERVIEW OF THE youre 


An overall schematic representation of the experimental 
apparatus is shown in Figure 2.1, and a photograph is shown 
In? Pei eee. The apparatus consists of a Pyrex-glass 
evaporator with the boiling test tube 1, a Pyrex-glass 
condenser 2, an R-114 liquid reservoir 3, an oil reservoir 
4, a graduated oil cylinder 5, a vacuum pump 6, a cooling 
section consisting of a 1/2-Ton R-l2 refrigeration unity 
tank filled with water-ethylene glycol peers 7, and a pump 
Ce The evaporator 1, condenser 2, oil-reservoir 4, 
graduated-oil cylinder 5 and R-114 reservoir 3, are placed 
in a Plexiglas chamber 9, where the temperature can be main- 
tained at’ lower “Vallee (up tcmeC mae relative to ambient 
temperature. This chamber also served as a safety barrier in 
case of an inadvertent overpressurization of the glass 
apparatus. 

The vapor rising from the evaporator 1s condensed in the 
condenser 2 and the liquid is then fed back into the evapo- 
rator through vals Wee In this manner, the specific 
enthalpy of condensation ("latent heat" ) of R-114 vapor is 
removed in the condenser. The evaporator of the R-12 
refrigeration system is made of a copper coil, which (ge 
immersed in the water-ethylene glycol tank in order to main- 
tain the temperature of the mixture at about -17 °C. An 
8-GPM, turbine-type pump 8 transports the water-ethylene 
glycol mixture thzeoenen the control valve VC to the 
condenser. 

The oil can be added into the evaporator 1 uSing the 


graduated cylinder 5 by opening valves Vl and V4. The 
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epraduated cylinder is supplied with oil by reservoir 4. The 
valves V5, V6 and V/ are provided for emptying the evapo- 
maeOr prior tomdismantling 1£, for example, diomechanging the 
evaporator tube. Also, fresh R-114 can be added to the 
reservoir while valve V5 is closed and V6, V/ are open and 
Peeproviding R-114 vapor through V8 and V1O. 

The test tube was placed horizontally in the evaporator 
after a 240 VAC electrical cartridge heater of 203.2 mm (8 
in.) in length had been inserted in its middle section. The 
whole system was connected to the vacuum pump through V8 and 
VlO in order to remove noncondensible gases. The relief 
valve was set to a gage pressure limit of 138 kN/m? (20 
psi), which is 50% less than the manufacturer-recommended 
working pressure limit of the Pyrex glass tees. All connec- 
tions were assembled with Swagelok fittings sealed with 
Teflon ferrules to ensure leak tightness. Also, all valves 
used in the system, except control valve VC and by-pass 
valve V9, were made of stainless steel. Copper tubes of 
9.5-mm (3/8 in.) diameter were used for all piping, except 


in the oil section, where 6.35-mm (1/4 in.) tubes were used. 


Eee SUOLLING TEST SECTION 


i, Evaporator 


The evaporator is a T-shaped container made of 


Corning Pyrex glass. Figure 2.3 shows a schematic of the 
evaporator with dimensions. The selection of Pyrex glass 
offers several advantages: it 1S corrosion-resistant and 


transparent and has a smooth interior surface (this mini- 
mized nucleate boiling at the inner surface of the 
container) compared to metals, and 1S more stable with 
temperature and pressure variations compared to ordinary 
glass. The operating pressure limit of the evaporator is 


guaranteed by the manufacturer up to a gage pressure of 200 


Zo 


KN/m-? “(Gps Each end of the Pyrex-glass evaporator was 
fitted with a cast-iron flange and a gasket. A detailed 
sketch of the cast-irom flange 1s Showman? cu poe eer Two 
aluminum flanges, 210 mm in diameter and 12.7 mm in thick- 
ness, were bolted to the cast-iron flanges. Thus, all 
fittings were connected through the aluminum flanges toueae 


Pyrex-ellass eCvaponaron. 
22. Lest Seciage 


A schematic drawing of the test section is shown in 
Figure 2.4. The boiling test tube was held in place by two 
Teflon bushings, which were attached to the aluminum flanges 
at both ends of the T-shaped evaporator. Four studs were 
used in order to attach the Teflon bushing on the aluminum 
flange. The Teflon bushing and both ends of the test tube 
were sealed by means of two rubber O-rings in order to 
achieve leak tightness. An additional O-ring was placed 


between the aluminum flange and the Teflon plug at each end. 


3. Boiling Tube 


A smooth, hard-copper tube, 15.9 mm (5/8 ieee 
outer diameter, 12.7 mm (1/2 in.) in inside diameétereand 
431.8 mm Cl? in.) \rmekene tine was used to provide baseline 
data for the nucleate pool boiling heat-transfer coefficient 
of pure R-114. The heater was a  1000-Watt 240-Volt 
stainless-steel cartridge, 6.35 mm (1/4 in.) in outer diam- 
eter ande203 .2 anne (Caen in length. The cartridge heater 
was inserted into a copper sleeve, which was 6.35 mm (1/4 
in.) in inside diameter, 12.7 mm (1/2 in.) in outside diam- 
eter and 203.2 mm (8 in.) in length. In order to provides 
Unit OF sea be bliss. the cartridge heater and the copper 
Sleeve were Soldered together. Pulido [Ref. 19] showed that 
(using a similar type of heater-element construction) the 


heater element provides a uniform heat flux along the 
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circumferential direction. The heater element used during 
Pas thesis was assumed to have similar characteristics. 
The cartridge heater and the copper sleeve were then 
inserted as a unit into the middle portion of the test tube 
Memme a tight mechanical fit. Gloteewr-wsiitimk-fit process 
1s essential to minimize the thermal contact resistance. 
However, owing to the time constraint in this thesis, the 
Shrink-fit process was delayed and a slide fit with a clear- 
ance of less than 0.01 mm was used. Thus, the resulting tube 
was used mainly for the purpose of obtaining a successfully- 
Operating apparatus). The active boiling length of the test 
mage was 203.2 mm (8 in.) in the middle portion of the 
Deming tube. In order to compute the actual average heat 
flux in the heated portion of the tube, a suitable correc- 
tion was applied for the natural convection created at both 
ends. See Heat Flux Calculation section in Chapter 4 for 
details. 

To measure the wall surface temperature of the 
boiling tube, 8 thermocouples were inserted into 8 grooves 
which were machined on the outside of the copper sleeve. As 
Shown in Figure 2.5, these thermocouples were located at 
different axial and circumferential locations. Four thermo- 
Ssempres (1,2,5,6) were placed with 90-degree separation 
around the copper sleeve to measure the temperature distri- 
bution at the mid cross section of the active length of 
bemlane tube. The thermocouple'’s wire was glued down at 
several places along the thermocouple's channel using Epoxy. 
The longitudinal temperature distribution on the active 
boiling section was also measured by the symmetrically- 
located 4 thermocouples (3,4,/7,8). The exact locations of 
the thermocouples and dimensions of the thermocouple grooves 
are given in Figure 2.5. All thermocouple grooves were 
axially machined from the location of the thermocouple hot 


junctions to the nearest end of the copper sleeve. 


Chel. 


C. CONDENSER. SECT IGy 


The condenser was also a T-shaped container made of 
Corning Pyrex glass. It was identical to the evaporator. The 
position of the condenser ‘can be seen in Figure 2.1. R-114 
was condensed on a vertical copper coil, which was inserted 
in the Pyrex-glass condenser. The copper coil was fabricated 
into a 76.2 mm (3 in.) diameter coil using 9.5 mm (37 Sie 
copper tube. The active condensation length was estimated to 
be 4.5 m (15 ft). | 

The top portion of the condenser waS connected to a 
portable, mechanical vacuum pump to remove noncondensible 
gases from the apparatus. The bottom of the condenser was 
also ponneeeed to the evaporator via valve V5 in order to 
return the condensed R-114 liquid to the evaporator. The 
cooling liquid, i1.e., water-ethylene glycol mixture, entered 
the top portion of the condenser, throtigh the copper @equm 
and left the condenser from the bottom to return to the 
water-ethylene glycol tank. The condenser was placed verti- 
cally and connected to the vapor outlet of the evaporator 
using L-shaped aliminum tube, 50.8 mm (2 in.) in diameters 
The maximum vapor velocity of R-114 vapor through the 
aluminum tube was found to be about 0.5 m/Sec. A Bourdog 
gage with a range of absolute vacuum to a gage pressure of 
1030 kKN/m* (150 psi) and a relief valve which was set to 138 
kN/m* (20 psi) were placed on the L-shaped aluminum tube. 


D. Ol ADeING SiGe 


To study the boiling performance of R-114-0il mixtures, 
a cylindrical aluminum reservoir, 152.4 mm (6.1n) “in doa 
eter and 152.4 em (6 in.) an height, and a glass "ou 
cylinder were installed above the evaporator. The relative 
positions of the oil reservoir andthe oil cylinder are 


SHOWN li iale epee The oil cylinder was 355 mm in lenge 
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and had a diameter of 25.4 mn. mits Cylinder was specially 
ordered to achieve a resolution of 0.5 ml. The oil cylinder 
was connected to the oil reservoir through valves V3 and V2. 
The addition of oil into the evaporator can be achieved 
through Vl by gravity after balancing the pressure of the 
o1l cylinder with that in the evaporator by opening valve 
V4. 


Eee CUOLING SECTION 


1. Water-Ethylene Glycol Mixture Tank 


In order to store the water-ethylene glycol mixture, 
a Special tank was manufactured. The total volume of the 
tank was 0.154 m* (0.48 m x 0.48 m x 0.66 m) and it was made 
of 12./-mm-thick Plexiglas sheet. All sides of the tank were 
glued together with methylene-chloride solution. The joints 
were held together with small screws for extra strength. The 
low thermal conductivity of Plexiglas was especially suited 
to minimize heat transfer (from room to water-ethylene 
glycol mixture) through the tank walls. The tank was placed 
on the floor and all sides were insulated with 22 mm (7/8 
mm thick rubber insulation sheets. The cooling mixture 
Semeained 49 tliters (13 egal.) of ethylene glycol and 94 
Mees Of (25 gal.) distilled water. The freezing point of 


mabeceemaxture was about - 25 °C. 


Yee R-12 Refrigeration Plant 


eee Ome h-Wizeret ri geration plant was installed to 
cool the water-ethylene glycol mixture. Figure 2.6 shows a 
mem@ematic OL the R-l2 refrigeration plant. It consists of a 
compact-type air-cooled condenser, a compressor, a receiver, 
a filter-drier unit, a pressure regulator, a pressure- 
control switch and a thermostatic expansion valve. The 


evaporator of the R-12 refrigeration plant was constructed 


des 


using a 9.5 mm (3/8 in.) copper tube, which was immersed saa 
the water-ethylene glycol tank. The temperature of the 
water-ethylene glycol mixture 1S controlled by both a ther- 
mostatic expansion valve and a pressure control switch. The 
R-12 refrigeration plant was adjusted to keep the tempera- 


ture of the cooling Ttqutd ae aroute ) eee 
3. Pump and Control Valve 


An 8 GPM, 115 VAC Burks turbine-type, positives 
displacement pump was installed on the floor and the 
Suction side of the pump, 25.4 mm (l in.) in diameter, was 
directly coupled to the water-ethylene glycol tank. Cooling 
liquid was pumped from the tank to the condenser through the 
control valve VC. Also, a by-pass valve V9 was placed before 
the control valve VC on the discharge line. The use of the 
by-pass line served two important purposes: 1, it avoided 
overloading of pump 8 in the event value VC is completely 
closed, and 2, it provided proper mixing for thew wage 
Stream returning from the condenser 2. See Figure 2.1 for 
the positions of control valve VC and by-pass valve V9. The 
by-pass valve was adjusted and set permanently so that a 
sufficient range of flow rates can be achieved through valve 


VC to satisfy proper cooling of the condenser. 


Pe Rooke Re oe Ole 


An aluminum cylindrical reservoir.) 22cm (2a in 
diameter and 254 mm (10 in.) in height, was placed verti- 
cally between the evaporator and condenser in order to store 
R-114 as a liquid. The liquid level of the R-114 can be 
observed by means of a sight hose attached on the reservoir 
with proper fittings. R-114 reservoir was connected to the 
vapor line through valve V/ and to the liquid line throuaa 


valve V6. See Figure 2.1 for arrangement of the reservoir. 
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G. CHAMBER 


Po alunimunetrame (1.0/7 mx O.51 m x 0.61 m) was 
constructed to locate all the parts of the apparatus, except 
the cooling section. All four vertical sides of the frame 
were covered with 12.7 mm (1/2 in.) thick Plexiglas sheets 
and both left and right sides were provided with hinges to 
enable easy access to the components of the apparatus. 
Aluminum and plywood plates were used to cover the bottom 
and top sides of the frame, respectively. The valve bodies 
of Vl through V8 were placed inside of the front Plexiglas 
Sheet with the valve stems penetrating the sheet. Thus, the 
valve handles were accessible from outside of the Plexiglas 
BOX . The whole frame was placed above the water-ethylene 
glycol tank with aluminum support so that the system was 
weway COMpact. 

One of the main advantages of this chamber is that the 
temperature Surrounding the evaporator can be reduced rela- 
tive to the ambient temperature. Also, in case of emergency, 
the thick Plexiglas chamber would provide a safety barrier 


to personnel and equipment. 


H. INSTRUMENTATION 
mee, Power Measurement 


A 240 Volt AC source was used as the power supply, 
and it was adjusted by a variac in the range of 0-260 Volt 
and 0-8 Ampere according to the desired heat flux at the 
surface of the boiling tube. Power input to the boiling tube 
was measured with an AC current sensor and an AC-DC true 
R.M.S converter. The AC current sensor was connected to the 
input line of the heater in series and the AC-DC true R.M.S 
converter was connected in parallel. Figure 2./ shows a 


Schematic representation of the power-measurement devices. 


or 


Both the AC current Sensor and the AC-DC (~Gmue Rei 
converter were connected to the data acquisition/reduction 


system. 
Pate Temperature Measurement 


Various temperatures were monitored throughout the 
system to include: 


Boiling tube wall (8 thermocouples) 


2. Liquid temperature (one thermocouple) 
3. Vapor temperature (two thermocouples), and 
4. Water-ethylene glycol mixture temperature (one ther- 


mocouple) 


The locations of the wall thermocouples in the sleeve of the 


boiling tube are shown in Figure 2.5. The liquid and vapor 
thermocouples were inserted LEG hie two specially- 
manufactured thermocouple wells. Figure 2.8 shows a sche- 


matic of these thermocouple wells. While the stainless-steel 
portion minimizes (owing to low thermal conductivity)" e11rome 
resulting from the axial conduction of heat from _ the 
Surrounding, the copper tip helps minimize the temperature 
drop from the area being measured to the thermocouple loca- 
tion (owing to the high thermal conductivity, Of scoppen es 

All the temperature measurementS were accomplished 
by 0.245 mm (30 gage) copper-constantan thermocouples. Each 
thermocouple measurement was read directly by a 
Hewlett-Packard 3497A data acquisition system, which was 
controlled by a Hewlett-Packard 9826 computer. Each thermo- 
couple was scanned for 0.8 seconds and twenty readings were 
averaged to obtain a more accurate measurement. 

A total of five thermocouples were calibrated. Two 
thermocouples were made from the beginning of a_ spool of 
copper-constantan wire; one was made from the mid portion 
and two were from the end. It was assumed that the proper- 


ties of the copper and constantan do not change along a 


a2 


peeveme SECtioOn of Wire for any given spool. All thermocouples 


were calibrated by the method described in Appendix D. 
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Figure 25 Schematic of the Boiling Test Apparatus. 
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a) Corning Pyrex Glass Evaporator (Dox D = 40 2xaaanme 
L = 178 mm, Ly,= 127 mn) 
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b) Cast Iron Flange and Gasket (d,= 190mm, dg= 210mm, 
L,= 14mm, A = ai 


Figure 2am schematic of the Pyrex Glass Evaporator. 
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III. DATA ACQUISITION / REDUCTION 


A. DATA ACQUISITION AND STORAGE 


A Hewlett-Packard 3497A automatic data acquisition/ 
control unit was used to read temperatures from the thermo- 
couples and to read current and voltage values of the 


cartridge heater from the AC current sensor and the AC-DC 


true RMS converter, respectively. A Hewlett-Packard 9826A 
computer was used to control the data acquisition/ conterem 
unit and to analyze and store data. Figure 3.1 showsmya 
photograph of the data acquisition/econtrey “unas and 


COmpuueer. 

Information was entered through the keyboard to prompt 
the data acquisition/control unit to automatically scan eaem 
channel. Channel assignments are listed in Table I. These 
raw data were immediately processed and a hard-copy printout 
was provided. Also, these data were transferred to a 
computer disk under a user-specified file mame to keep a 
permanent record. The ability to store raw data directly 
enabled these data to be reduced at any time and allowed 


flexibility for changes to data-reduction software. 


B. DATA REDUCTION 


Following data acquisition for each data point, resimilme 
were computed according to the stepwise procedure outlined 
in the next section, and then printed on a Hewlett-Packard 
26/71G thermal printer. Heat flux versus wall superheat (the 
temperature of the boiling surface minus the liquid tempera- 
ture) were also stored in a user-specified plot data file 
for subsequent plotting using the subprogram PLOT. A 
Hewlett-Packard 7470A plotter was also interfaced with the 


COmeouie ease 
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Sue WT SteeATA-COLLECTION AND SOLUTION PROCEDURE 


He 


HO: . 
cl. 


eZ. 


iS. 


14. 


Select tube type (all dimensions of the boiling test 
tube are included). 

Set desired heat flux and saturation temperature of 
Pec womenien Liduid. Waite Goresteady-state conditions 
(See Chapter 4 section B for details). 

Scan all channels listed in Table I (thermocouples, 
AC current sensor and AC-DC true RMS converter read- 
aie. S je 

Save these raw data in a user-specified file. 

Convert these raw data readings to corresponding 
units (temperature in °€, current in Ampere, voltage 
id Volt ) 

Compute the heat-transfer rate from the cartridge 
heater (See Chapter 4, Section C for details). 
Compute the average wall temperature of the boiling 
test tube and calculate the wall superheat of this 
data set (See Chapter 4, Section C for details). 
Compute the physical properties of R-114 using given 
correlations at film temperature (See Appendix H for 
details). 

Compute the natural-convection heat-transfer coeffi- 
cient of R-114 from non-boiling ends of the test 
tube. 

Compute heat losses from non-boiling ends. 

Calculate the heat flux from boiling test tube to the 
boiling liquid. 

Calculate boiling heat-transfer coefficient of the 
R-114 from test tube. 

store the heat flux versus wall superheat values for 
each data set in a user-specified plot file. 


Use the subprogram PLOT and plot the data run. 
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Appendix E shows a listing of computer program DRPRI1 and 
subprogram PLOT and Appendix F shows an example of represen- 


tativemedaka arnum. 


TABLE I 
HP 3497A Channel Assignments 


Channel Assignment 

2 eS 2 Boiling tube wall temperature 
35 Liquid temperature 

34 - 35 | Vapor temperature 

56 Temperature of cooling liquid 
62 AC-DC true RSM converter 

63 AC current sensor 
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Photograph of Data Acquisition/Reduction System. 
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IV. BABPERIMENTAL SERUCe DURE 


A SP REA AT TON 


1. Pressure Testyotmeme Apparatusmymen Air 








Upon assembling all the components, includingwehe 
boiling test tube, the experimental apparatus was pressur- 
ized With air up to TOG (15 psi A soap-bubble test 


was first carried out to locate any leaks from each compo- 
nent of the apparatus. All detected leaks were succesfully 


and systematically fixed at the end of each presSure test. 


2. Pressure Test of the Apparatus with R-114 








A second pressure test was performed using R-114 
vapor. The apparatus were charged with R-114 vapor in the 
following manner. The apparatus was evacuated to 25 in. Hg. 
by a portable mechanical vacuum pump through valves V1ll and 
V8 (See Figure 2.1 for the configuration of valves). Upem 
c leis i nee al em ee a vapor outlet of the R-114 supply 
cylinder was connected to the apparatus by means of valve 
VlQ. Thus, the pressure of the apparatus was increased up to 
the saturation pressure of R-114 at that ambient temperature 
(ambient temperature was 21 °C and the corresponding satura- 
tlon pressure (gage) was 86 kN/m* (12.5 psi)). An Automatic 
Halogen Leak Detector, TIF 5000, was used to detect R-114 
leakage from the apparatus. The sensitivity of this detector 
1 Sly 2 os per year, 3 PPM max. A few small leaks were 


observed and successfully isolated. 


3. Charging the Apparatus with R-114 


Upon completing the pressure tests, the evaporator 


of the apparatus was filled with R-114 liguid up to a mariges 
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position (20 mm above the boiling test tube) with the 


following procedure (See also Figure 2.1). 


Mis 


Be 


The temperature of the water-ethylene glycol tank was 
reduced to about -17 °C by the R-12 refrigeration 
else 


The apparatus was -evacuated to about 29 in Hg vacuum 


.by the vacuum pump (the valves Vll and V8 were open 


and V1O was closed). 

Voce los tae VY 11, the vapor outlet of the R-114 
Supply cylinder was connected to VIO. 

The cooling section pump 8, delivered the cooling 
liquid frome the water-ethylene Bimeoeumcamie Fo the 
condenser through the control valve VC. 

The R-114 vapor condensed on the condenser coils and 
the liquid was collected in the evaporator by 
eravity, while V5 was open and V6 was closed. 

Pump 8 was then stopped when the desired liquid level 
of the evaporator was achieved. 

The pressure of the apparatus was then allowed to 
increase up to the saturation pressure corresponding 


to the ambient temperature. 


NORMAL OPERATION 


The following procedure was established to obtain the 


heat-transfer coefficient of R-114 from the smooth copper 


test tube: 
1. The test tube was immersed in the pool (20 mm below 
the liquid level). 
meee tne R-1l2 refrigeration unit was operated 24 hours in 
advance in order to reduce the temperature of the 
water-ethylene glycol tank to about -17 °C. 
rl neomGdomandeGqud@scitilon,control unit, computer, power 


supply and cooling section pump were switched on. 
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OR 


Any accumulated moncondensible gases were evacuated 
by the portable mechanical vacuum pump through valve 
V8 and V1ll (see Figure 2.1). 

The iteractive computer program DRPR1]1 was operated 
and the desired heat flux and saturation temperature 
were given to the program as reference values. 

The cartridge-heater voltage was then adjusted by 
variac in order to maintain the desired heat flux. 
The desired and actual heat fluxes were compared 
continiously by the computer program until they 
agreed to within 2 percent. 

The actual versus desired saturation temperature of 
the liquid R-114 was also monitored by computer 
program. The amount of ecoot nie aura. which was 
being circulated through the copper coil in the 
condenser, was regulated by control valve VC in order 
to obtain nearly constant Saturation temperature 
(40.2 °C) at a leevemeie tele 

For all the consecutive Settings during both 
increasing and decreasing heat fluxes, the boiling 
was allowed to stabilize for five minutes at each 
power setting and:‘raw data (thermocouple readings, AC 
current sensor readings and AC-DC true RMS converter 
readings) were recorded in a user-specified file. 

At a given power setting anda saturation tempera- 


ture, the following processed data were recorded as a 


printout: wall temperatures of the boiling section, 
liquid bulk temperature, vapor temperature, net heat 
£ dies temperature of the water-ethylene glycol 


mixture, wall superheat and heat-transfer coefficient 
of the R-114. 
For each data set, the above procedure beginning with 


Step 5 was repeated. 
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ii. For each data series, seven different heat fluxes (1, 
2, 5, 10, 20, 50 and 95 kW/m?) were selected. 


C. HEAT-FLUX CALCULATION 


Mmacordinee Geunthne deSeription Some DoOLiing tube 
construction, the cartridge heater, which was inserted into 
the boiling tube is the heat source. The variations in heat 


flux are made through different voltage settings of the 


Varlac. 
e = VI oe) 
where: 
ee oe ea ee eo et Lome tame ycartridge 
heater (W) 
Ve mevGleare across the cartridge heater 


element (volts) 


P= current through the heater element (amps) 


Ms Clrcul1t 15 connected in series with an AC current 
Sensor and in parallel with an AC-DC true RMS converter. 
Each of these two sensors produces a DC output in the range 
wee -10 VV, which iS automatically read by the data 
@eemrsition/control unit. The calibration equations of these 


Sensors give the actual values as follows: 


V = 25V- (2) 


ce), i (3) 
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where: 


V_= voltage reading by AC-DC true RMS 


converter (volts) 


I. = current reading by AC current 


sensor (volts) 


The total length of the boiling test tube was 355 ¢iimm 
while the active boiling length; i.e., the length of the 
cartridge heater, was 203.2 mm. The geometry of the boiling 


test tube is shown schemativeally wGneneouno ee 


EPOXY BOILING COPPER 
Pen WE Sea ve 





Figure 4.1 Geometry of the Boiling Test Tube. 


Natural-convection heat transfer was assumed to occur at 
both ends of the test tube (the distance LL). This paraa 
the tube was considered to be a Straight fin of uniform 


cross section. For this cylindrical geometry: 


Al = 7 (Dg - Dt)/4 (4) 


C 


DO 


mwmere : 


A. = pL (2) 


cross-sectional area of tube (m7?) 
tube outside diameter (m?’) 
tube inside diameter (m’) 


tube outside surface area of 
nonboiling section (m7?) 

tube outside wall perimeter (m) 
i 


O 


non-boiling length of the test tube (m) 


It was assumed that the temperature at the base of the 


Stieaight fin 


was equal to the average wall temperature of 


the active boiling section. The average outer wall tempera- 


ture of the active boiling section was calculated using the 


radial conduction equation from the thermocouple position to 


Smmerstirface of the tube. 


wimere : 


8 
ae () Ta) /8 (6) 
Twot Taye Q( IN(/ DI/(2TL KE) ) (7) 


Sue ce 


temperature of the thermocouple 


el 


Se ce 


BReyerch calico (Gel): 


Tavg average wall temperature at the 
thermocouple locations (°C) 

D, = diameter sae) the posiritem Or Ehe 
thermocouple (m) 

D = outer diameter of Ehe Beriing tube Gay 

L = active boiling tube length (m) 

k, = thermal conductivitiy of the copper (W/mK) 


outer wall temperature of the boiling 
test) Cuber(-c) 
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T, = temperature at the base of straight fie 


The temperature distribution of a fin of uniltlormyenre 2 


ional area is given by [Ref. 20] with the follows 


assumptions: 


ile One-dimensional conduction in the x direct ion. 
2. Steady-state condition. 
3.  @onstane therma  conauweriv 1 awe 
G. Negligible radiation from the surface of the fin. 
5. No heat generation in the fin. 
6. Uniform convective heat-transfer coefficient over tame 
Ei ieestiteat ace. 
d*6@ / -dx= ise m=an0 (8) 
where: 
Q(x) = T(x) - Teat (9) 


DZ 


and 


io) ae py ko Ac (10) 


where: 
a eo mee om wemperadaure on fluid eG.) 
hEEEEECOnecetaive heat-transter 
coefficient (W/m?.K) 
Now, assuming negligible heat loss from the thin tube 


fee at x = L,,) 
(de / haa 0 Cine) 


peat sSO at the base: 


I 


Qs eee ee ete Lo. 
where: 


Qp= heat-transfer rate from the base (W) 


Eeine the above boundary conditions (Eq. 11, 12), the 
temperature distribution along the straight fin can be 


expressed as: 


(0/O,) = ( Cosh(m(L,, - Seay ees) Cosh(mL_ ) ) cise) 
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and the fin heat-—t vamoremmneaecwa .. 


Qe = Qpzaf/h p kA. OpTanh(mL, ) exs 


Or = heat-transfer rate through one non-boiling 


length of the test tube (W) 


where: 


The average difference between the wall temperature and 


the liquid temperature may now be determined from. the 


following equation: 


by 
(98/98,) = ([(cosn(mt,- x) )/Cosh(mL, ) Meee LS @bs) 
6 = fo mLy) Tanh(mL,) 
6) = iyo ae 
where: 
aes Lat average diitference betvcensouten 


wall temperature and liquid saturation 


temperature (K) 


A free convection correlation stated by Churchill and 


Chu [Ref. 20] was applied for the average Nusselt number on 
a horizontal cylinder: 


mic 
0.387 Ray ) 2 a 
Nu = 10.60+ - = 
oS | C1 2 (OSS 07 ete 


107° < Ray < 10° 
O 
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The average Nusselt number ls: 





Nu, = k 
O 
where: 
k = thermal conductivity of R-114 (W/m.K) 
Solving for h: 
Pec (18) 
mm i | D % 
h = — (0.60+ sie 
o | [1 + (0.559/pr) 774°) 8/2! | 
Ce) 
Pr = = 
a 
a 2 2.0 
= — IBMT Taq) Po ee 
D 7 v a 
O \ 
Oeeeacl 
where: 
Je = Weskercle} iolviulersne 
Vv = kinematic viscosity (m?*/s) 
® = thermal diffusivity (m?/s) 
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Rano = average Rayleigh number 
g = gravitational acceleration (m/s?) 
P = volumetric thermal expansion 


coefficient (1/K) 


© = density of R-114 liquid (kg/m?) 


Now, substitution of Equations 10, 15, 20 into Eqiuaiaaae 


18. tesa hese: 








— 72 
gsD-e Tani 2 Lu | ae 
oOo B Uk AG 
Z 
ieee 2 

: oo (220) 

h = —<0.60 +0.387 = 

Oe [ie NOMNeCo Per 
Equation 22 is solved for h by an iterative technmeme 
within a range of precision of Q.001. Knowing the value of 


the natural convection heat-transfer coefficient along the 
non-boiling ends, the total heat-loss rate is calculated 
from Equation 14, and the heat-transfer rate through the 
boiling section is obtained by subtracting the total heat- 


loss rate Erom the total hear-trans renner 


° 5 (23) 


Q = Og = 2xQ, (24) 
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Haemineat Clix from the boiling surface is: 


S| Ss Oe (22) 


oe 





V. RESULTS AND DistUas ron 


A. OUTLINE, OF THE DATA 


Using the procedure outlined in Chapter 3, nine data 
runs were completed, primarily to de-bug the experimental 
apparatus for its successful operation. For all data xunce 
the same smooth copper tube, described in Chapter 2, was 
used as the boiling surface. Each data run consisted of 
seven different heat fluxes (1000, 2000, 5000, 10000, 20000, 
50000 and 95000 W/m*) with a specified saturation tempera- 


ture S02) slemor 20RCc = All the raw data were Stored Gime 
user-specified file, named "“WHxx”, where "W" indicates 
Wieland tube, "H" indicates hard copper and "xx" indicates 
the data run number. Also, the heat-transfer coefficient 


versus heat flux data were stored ina plot file named 


tf 


OP ex A summary of these data runs is piven in Tab leu 


B. LONGITUDINAL AND CIRCUMFERENTIAL TEMPERATURE VARIATIONS 


During all data runs, especially at high heat fluxes 
(greater than 10 kW/m?), large temperature variations were 
observed along the active boiling length of the test tube. 
Figure 5.1 shows the measured wall temperatures at each of 
the thermocouple locations, while Figure 5.2 shows’ the 
temperature variations in the circumferential direction (at 
the axial mid point) at two different heat fluxes (2 and 92 
kW/m?*). Despite the uniform heat flux provided by the 
heater, the temperature measurements Showed considerable 
variations in both axial and circumferential directions 
The reason for these temperature variations may be explained 
as follows: As explained in Chapter 2, during the manufac- 


turing of the boiling test tube, the copper sleeve was 
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/Mioerted into the boiling test tube without any shrink-fit 
process. The slip-fit process (with a diameter clearance of 
about 0.01 mm) used inthis study was believed to have 
caused considerable thermal contact resistance at the inter- 
face between the boiling tube and the copper sleeve. Also, 
it is quite possible that this contact resistance was not 
Moeroarm alone both the axial and circumferential directions. 
The circumferential temperature distribution may also be 
Seer iputed to a2 circumferential variation in boiling heat- 
transfer coefficient. Ditapnca tives higin neat flux boiling 
regime, it was observed that the rate of bubble formation at 
the top of the horizontal boiling tube was lower than at the 
rest of the boiling tube; 1.e., the boiling heat-transfer 
coefficient at top of the boiling tube was relatively lower 
than the other locations. This poorer performance at the top 
of the tube may be attributed to the fact that this portion 
was mot receiving sufficient amount of fresh liquid to 
replace the liquid that evaporated. [new ob cieenectlLonm Of thie 
liquid flow could be caused by the growing vapor blanket 
from the bottom to the top of the tube. Based on this obser- 
vation, it is expected that the temperature at top of the 
boiling tube could be greater than the other locations. As 
Shown in Figure 5.2, location 5, which represents the top of 
the horizontal boiling tube, measured the highest tempera- 
ture difference as expected from the above-mentioned obser- 
Weenie LOT) . 

However, positions 1 and 2 showed a sizable discrepancy 
at high heat flux, and this was due most likely to the non- 


Bumenormm Contact resistance. Also, a comparison of thermocou- 


ples 2, 7 and 8 show a similar discrepancy. It 1s clear 
therefore, that the shrink-fit process must be used in the 
future. 


a 


C.-PLOT ANALYSIS OQRBRMGCLEATES SOL ee 2G ee 


Figure 5.3 Shows a typical nucleate pool-boiling 
performance curve of the smooth copper tube in R-114. The 
observed behavior of this process is analyzed by studying 
the different regions Gf 7Ehe boiling veurm Ee. 

From point A to pe ime. . a continuous increase in wall 
superheat (Ge, = Te) 1s observed when heat flux is increased. 
Also, during the experimental runs, no bubbles were observed 
along the boiling tube in this region of the curve. This 
region corresponds to the natural-convection process. 

ELOm pOLDt bY tO pone. a reduction in wall superheat 
is observed, while the heat flux continuously increased. 
While the middle portion of boiling tube showed bubbles, the 


rest of the boiling tube showed no bubbles during the exper- 


iment in this region. This region is known as7~ the mixed 
boiling region, where transition from natural convection to 
nucleate pool boiling heat transfer takes place. The wall 


Superheat continued to decrease until numerous nucleation 
Sites became active (point C). 

After point C, the wall superheat starts to increase 
with increasing heat flux as shown by region C to D. A very 
high density of bubble formation was observed along the 
active boiling section of the test tube (only a few bubbles 
were observed on the nonboiling sections of the test tube). 
See Figure 5.4. This observation confirmed the validity of 
the natural convection model assumed for data reduction. 
However, the R-114 liquid returning to the evaporator from 
the condenser caused a mild turbulance in the entrance 
region of this return line in the evaporator. Inherently, 
this affects the natural-convective heat-transfer perform- 
ance of nonboiling ends of the test tube. 

When heat flux is gradually decreased, the curve follows 


a different path (from D to E) as shown im Figure 5.37 ) Gime 
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is due to the existing, stable nucleation sites remaining 
active for a wide range of heat fluxes. Similar analysis can 
also explain Figure 5.5, which represents the same data run 
Puimom a ditterent basis {(heat-transfer coefficient versus 
reat flux). 


Pee OP RODUCIBILITY TEST OF THE APPARATUS 


In order to test the reproducibility of the experimental 
apparatus, two runs were performed on two different days 
(data run number 5 and 7 in Table II) at the same saturation 
temperature. Figure 5.6 and Figure 5./7 show the comparison 
of these two data runs. It can be seen from these figures, 
that there is very good | agreement between these two 
different runs. This agreement shows the ability of the 
apparatus to reproduce data runs revealing successful 


Seeration. 


PemeelOLLING PERFORMANCE OF SMOOTH COPPER TUBE IN R-114 


In order to test the validity of the data taken from the 
present experimental apparatus (despite the presence of 
errors owing to the thermal contact resistance mentioned 
earlier), an attempt was made to compare the present data 
with data found in the literature. For this purpose, data 
memumber 5 (10 °C saturation temperature with decreasing . 
heat-flux rate) was compared with two sources in the litera- 


meme as discussed below. 


i Comparison with Chongrungreong-Sauer Correlation 


A correlation proposed by Chongrungreong and Sauer 
[Ref. 7] for the nucleate boiling performance of refrige- 
mames and retrigerant-oil mixtures 1S compared with run 
number 5. This correlation is based on a dimensional anal- 
feeseand Gata from various sources: [Ret. 7), (Ref. 14]; 
Peres loSigmukreft, 22) and [Ref. 23]. 
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Aaa? S69 


PA hig 





ah 395 (26 | 


x [pyres | D ~( 444 a 1 579 | 
10.01568 aan 





where: 
h = heat-transfer coefficient (W/m?.K) 


Q = heat rate (W) 


A = surface area (m7?) 
D = characteristic dimension of heated surface (im) 
WH, = viscosity of saturated liquid (@r/mee 


oe Tat ene heat of vaporization (W.s/g) 
C.= specific heat of MiquitanGy er. 
k,= thermal conductivity of liquid (W/m.K) 
P = pressure (atmospheres) 
Q. = volume fraction of pure refrigerant 
OS eet Ut leche c rnin (g/cm?) 


Op. = vapor density 2/cem | 


Eq'n (26) agrees very well with theirdata [Ref. 7] 
for R-1ll as well as with Stephan s data [Ref. 22) forsoueee 
free refrigerants, R-ll, R-12, R-13, R-21L, R-22, R-l US 


R-114 with errors smaller than 16 percent. Figure 5.8 pres- 





ents a comparison of the current experimental data with eqin 
(26). As can be seen from this figure, the current data are 
in excellent agreement with this predictive equation. (The 
analysis given in Appendix I shows that the uncertainty in 
the boiling heat-transfer coefficient for this experimental 


investigation is about 20%.) This agreement was somewhat 
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unexpected since the present data includes some kind of 
contact resistance in it. If this contact resistance were 
known, and were removed from the data, a higher boiling 


@@elficient would result. 
Eee COMparisem with Data of Henrici 


The nucleate boiling heat-transfer coefficient of 


R-114 and R-114-0il mixtures from a horizontal, smooth 
copper tube (30 mm in outer diameter) has been experimen- 
tally investigated by Henrici [Ref. 10]. Since Henrici’'s 


data were taken from a tube with a larger diameter than the 


current boiling tube diameter, his data were revised _ to 
include this difference. For this purpose, a semi-emprical 
correlation to show the effect of tube diameter, as devel- 


oped by Cornwell et al. [Ref. 21], was used: 


itl =e b onamaas Ca 


where: 
Re = (q D) / a and Demeenchey / k 


C =iSOefor reifigerant 


Using this relationship, Henrici's data were revised 
moa «=€6diaAmeter of 15.875 mm (5/8 in.) used in this 


experiment. 


ee om onum mein 30emm) (30/15.875)*?? 


h(d=15.875 mm) = 1.24 h(D=30 mm) 
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A line representing the revised Henrici’s data is 
plotted in Figure 5.9 together with the unrevised data line. 
It can be seen that the boiling heat-transfer coefficient 
increases by about 25 percent when the diameter is decreased 
from 30 mm to 15.875 mm (50-percent reduction in diameter). 

Figure 5.10 provides Henrici Ss natural Convectmwen 
data and Henrici'’s revised boiling data in comparison with 
the present experimental data. As expected, the natural- 
convective heat-transfer coefficient is lower than the 
boiling heat-transfer coefficient. The boiling heat-transfer 


coefficient from Henrici's data is higher than the coeffi- 


client found, ait hiseeckea, The reason for this may be 
explained as follows. In Henrici's experimental investiga- 
tion, the measurement of the surface temperature has been 


done uSing thermocouples outside of the boiling surface. 


Unfortunately, he did not describe in detail how he measured 


the wall temperatures. If the thermocouples were not 
totally imbedded in the wall, they may produce incorrect 
wall temperature values. If a portion of a thermocouple is 
in direct contact with the liquid, this thermocouple waa 


measure a value lower than the actual wall temperature and 
such an error would decrease the wall superheat, thus 
resulting in a higher heat-transfer coefficient. Also, "a5 
the contact resistance effect 1S removed from the present 
data then the heat transfer coefficient would increase 


toward the Henrici data. 


BP. EPEPEC USO Sebo ouUR: 


The effect of pressure on the nucleate-boiling heat- 
transfer performance of R-114 was investigated for reduced 
pressures (2) Ee Ne “Of SOR 0Z 6a. 0.0391 and 0.0556 (data maa 
numbers 5, 8 and 9, respectively in Table II). Figure 5 7 


Shows the relative increase in boiling heat-transfer 
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coefficient with increase of pressure as expected. In order 
to calculate the pressure effect on the boiling performance, 


a simple equation developed by [Ref. 7] was used. 


hemo? (q)e°° (ONees p24 (28) 
where: 
h = heat-transfer coefficient (W/m’.K) 
q = heat flux (W/m?) 
Om=—/VOMunNeMiraetlen of pure refrigerant 
P = boiling pressure (atmospheres) 
Using this equation, the pressure effect on the boiling 


femaeeerransier coefficient can be written as follows: 


h(P,) = h(P,) (P, /P, )%?5 


According to the above equation, the effect of pressure 
from 127 kKN/m* to 87 kN/m? decreases the boiling heat- 
transfer coefficient by about ten percent, which agrees 


reasonably well with the current experimental data. 
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Figure 5.2 Circumferential Temperature Variation 
on the Boiling Tube. 
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VI. CONCLUSIONS 


An experimental apparatus to study the nucleate- 
boiling heat-transfer performance of R-114 from a 
Single horizontal tube has been designed, constructed 
and instrumented. 

An analysis scheme to establish the nucleate boiling 
heat-transfer coefficient of R-114 has been outlined, 
and data acquisition/reduction programs have been 
written and tested. 

The experimental apparatus was successfully operated 
in the saturation temperature range of 0 °C to 20m 
and in the heat flux range of 1 to 95 kW/m’. 

The test runs performed on different days on the 
apparatus showed repeatabilty within + 1 percent, 
revealing successful operation of the apparatus. 

The uncertainty in wall superhet (about 20%) drasti- 
cally affected the boiling heat-transfer coefficient. 
This can be attributed to the large temperature vari- 
ations measured along and around the boiling tube 
because of a non-uniform contact pressure between the 
copper sleeve and the boiling tube. 

The data taken on this experimental apparatus have 
been compared with the experimental correlation for 
refrigerant and refrigerant-oil mixtures developed by 
Chongrungreong and Sauer [Ref. 7]. The present data 
agreed to within 5 percent with this correlation. 

The current data were also compared with Henricigg 
(Ref. 10] experimental nucleate boiling data (from 
Smooth horizontal copper tube in R-114). The present 
data lie as much as 50 percent below Henrici's data. 


This disagreement may be attributable to the thermal 


Ve 


eem@eacth  westormanice present dm the bolling tube, or to 
the errors in Henrici's wall temperature 


measurements. 
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VII. RECOMMENDATIONS 


To prevent considerable temperature variations along 
and around the boiling tube, a $shrink-fit process 
must be applied in order to have more uniform contact 
pressure between the copper sleeve and the boiling 
tube. 

To determine the effect of contact resistance 
(between copper sleeve and boiling tube) on _ the 
nucleate-~boiling heat-transfer performance of a hori- 
zontal tube, a series of experimental investigations 
should be performed with different contact pressures. 
To maintain constant electrical power input into the 
boiling tube, a voltage regulator should be added to 
the input line. 

To keep the temperature of the water-ethylene glycol 
mixture at a constant value, the 1/2-ton Capacme 
R-12 refrigeration unit should be replaced with a 
IL Teteng ibhglabie . 

The positive-displacement pump should be replaced 
with a centrifugal-type pump, in order to deliver 
lower flow rates without reaching an electrical over- 
load limit. 

To check the saturation temperature of the boiling 
liquid corresponding to the measured saturation pres- 
Sure, amore accurate pressure transducer should be 
used. 

The active length of the cartridge heater should be 


carefully measured or obtained from the manufacturer. 
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APPENDIX A 
PHYSICAL PROPERTIES OF FREON FLUOROCARBON COMPOUNDS 


“FREON” 11 | “FREON 12 “FREON” 113 “FREOM” 114 
Enemies or muls CCIE CCIF, CCIF, CCIF, 
Molecular Weight 13737 120 92 187 38 170 93 


Boiling Point at 1 atm i 23 82 — 2979 4757 S77 
or 7487 — 21 62 117 63 38 78 

Freezing Point ac — 158 5 = 04 
— 252 —jl —137 


tr 
1120 214% 
4174 










Critical Temperature mG 
388 4 


Critical Pressure atm 435 ayy 
Ibs/sq In abs 639 5 495 
Critical volume ec/mol 247 325 
eo ft/ib 00289 0 0278 
Critical Density g/cc 0.554 0 558 0.576 
Ibs/eu ft 346 348 360 











Density, Liquid g/ce 
at 25°C (77°F) ibs/eu ft 
at 54.5°C (130°F) g/cc 

Ibs/cu ft 





Density, Sat'd Vapor g/l 5 86 633 

at Bosling Pont fhs/eu ft 0.367 0 395 
Specific Heat, Liquid 

(Heat Capacity) cal/(g) (°C) 0 7232 0218 0243 
Speciffe Heat, Vapor, 

at Const Pressure (1 atm) ~~ cai/(g) (°C) 38°C 60°C 

at 25°C (77°F) Btusciny CAYO @ oor 0.161 @ (149°) 0170 


at 25°C (77°F) Btu/(!b) (“F) 
Specific Heat Ratio of Vapor 


at 25°C and | atm Cp/Cy 





Heat of Vaporization cai/g 43 10 
at Boiling Point Btuy/ib 77.51 





Thermal Conductivity at 25°C (77°F) 
Btu/(hr) (ft) (°F) 
Liquid * 0 0506 0 0405 00434 
Vapor (1 atm) 0.0045] 0.00576 00014 
(Data from ASHRAE in most cases) (0 5 atm) 
Viscosity at 25°C (77°F) 
Liquid centipoisa 0 430 0.214 0.68 0 36 
Vapor (1 atm) centipoisa 0.0107 00123 0910 00112 
(Data from ASHRAE In most cases) (0 1 afm) 
9 


Surface Tension at 25°C (77°F) 173 12 
dynessem 
374 1.287 ] 354 1 288 
2.46 








Refractlya index of Uquid 
atz5 C (77°F) 


] 
Retative Dielectric Strength of Vapor al | atm 3.71 
and 25°C (77°F) (nitrogen = 1) 


39 (0 44 atm) 334 


Dielectric Constant 


Liquid 228 @ 29°C 24] @ 25°C 226 @ 
Vapor (1 atm)*s 1.0036 @ 24°C% 1.0043 @ : 
0017 
(Sat'n Pres) 





Solubility of Water in “Freon” 





Solubility of “Freon” In Watar 
at ] atm and 25°C (77°F) wt % 


at 25°C (77°F) wt % 
Sotubility Parameter ( ) 75 72 62 
Kauri-Butanol Vaiua (KB) 60 32 12 
Toxicity’ pom (v/v) 1000 1000 1000 1000 
Threshold Umit 

Vaiua (TLY) mg/m’ 5600 4950 7600 7000 
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APPENDIX C 


A | 


APPLICATIONS OF "FREON'' FLUOROCARBON COMPOUNDS 


, The tabfe below is intended to provide a general view of the range of applications and is not all inclusive. 


For specialized applications of more detail, please make specilic request. 


Solvents, Blowing Agents, 
Fire Extingquishants, 





applications, usually in “cascade” 


Ftuorocarbon Refrigerants Aerosol Propeltantsib!) Dielectric Fluids and Other Uses 
“Frean’ 23 Component of ‘Freon 503 

azeotropelal 
7 saan pe 


systenis 





Dielectric fluid 


Intermediate berween ‘Freon’ 13 Efficient fire extinguishant 

and “Freon” 22 lor medium to low (Haion' 1301) especially surted lor 

temperature applications Not automatic protection of materials 

extensively used. subject to water damage and of areas 
occupied by personnel 


“Freon” 22 Household and commercial 
refrigeration and ait conditioning 
anplicatians, Perinets use ol snvailer 
equininent. Component ol 


azeotropes!3). 













Orelectric thiuid, an economic 
replacernent for “Freon” 116 in 
Most dielectric applications. 


Accepted as a food propellant by the 
FOA, this material is weil suited for 
food aerosols and finds use in fat 
emulsion food whips Good foam 
stability with absence ol odor or taste. 


Used as an azeotrope component in 
“Freon 502!8). 














Blowing agent lor foained plastics 
applicationsid}. Qvelectric gas 


Most widely used high pressure 
propellant fur non lood use Blends 
with “Frean” 11 and “Freon” 114 
are widely used. 


Most widely used relrigerant in 
household, autontuntive and 
commercial relrigeration and air 
conditioning sysiemis Also as a 
component of azeotropes!a! and, 
in high purity lorem (“Freon” 
freezant} approved as a direct 
cuntact lreezing agent for foods. 







































Low nressure propellant, aiternative Blowing agent for loamed plastics. 
to “Freon” 11, having poorer 
sotitbility properties and less odor. 


Especially used in personal products 


a eat transter ttn 


“Freon” VI Widely used in centrifugal Most widely used low pressure Occasionally used as a solvent 


compressors for industrial and propellant for non-food use Ooes (Freon MF) Biowing agent for 
commercial are conditioning systems Not provide adequate plessure alone. foamed plastics! 
and tor industrial cooling of process $0 1s almost entrrely used in blends 
“Freon” 113 
{al A number of azealrapes I’ Freon’ 500, Freon 50? etc | are available tor refrigeration use Bulletins describing the camposttion, 
Properties and uses Of these mixtures ere available on request. 


in farge industrial process cooling 
and air conditioning systems using 
multi stage centrilugal conipressors. 
































water or brine Low viscostty and with “Freon” 2c), 
freezing point permit use as a 
fow temperature cooling liquid. 
























Solvent in some aerosol! lormulations, 
usually propelled with “Freon” 12 


fn commercial and industrial ave 
conditioning and process water or 
brine chilling using centrilugal 
Compressors, particulary in smalf 
tonnage applications. 


Extensively as a solvent (“Freon” TF} 
alone and sn special purpose 
formulations for a wide range of 
critical cleaning needs In cutting 
fluid formulations, Vaiclene 

dry cleaning, etc. 











Ibt Normally aerosol propellants are blended to qve the required vapor pressure and solubdily requirements 
Ic] “Freon” 11 Sis a siatiized grade trequentty used in formulations containing alcohols or waler 


(d) “Freon” t2 and “Freon 11 are widely usec as blowing agents for a range of foamed plastics in which [hey Drovide excellent cell structure 
and Ihe trapped Huorocarthon in closed cell loams siqnilicantly imoroves insulation and recovery properties Ot the loam 
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APPENDIX D 
THERMOCOUPLE CALIBRATION 


A. EQUREMENT Use. 


The equipment used in thermocouple calibration is shown 
in Figure D.1l, and a brief description of each component is 


given below: 


1. Thermocouple Wire 


Type-T (copper-constantan) Teflon-coated wire of 
0.254 mm (0.01 in.) in diameter was used for all thermocou- 


ples. 
2: MG@alibprats onesie 


A Thermos flask was used as the calibration bath. In 
order to maintain an isothermal temperature distribution, a 
motor-driven mixer was used. Observations during thes ca. 
bration procedure showed bath temperature fluctuations to be 
+ 0. CC 2B 


3. Thermocouple Readout 


A Hewlett Packard 3497A automatic data acquisition/ 
control system and a Hewlett Packard 9826 computer were used 


to read, analyze and record the calibration data. 
4. Reference Temperature 


A Hewlett Packard 2804A quartz thermometer was used 
to measure the bath temperature. This quartz thermometer had 
a resolution of 0.0001 K, while the manufacturer-guaranteed 


accuracy was better than + 0.03 K. 
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B. PREPARATION FOR CALIBRATION 


1. Thermocouple Preparation 


A total of five thermocouples were prepared for the 
eee pration run. Two thermocouples were made from _ the 
beginning of a spool of copper-constantan wire; two were 
made from the end; and one was made at the mid section using 
the following procedure: 

* The Teflon insulation was removed for a length of about 
4 mm from one end of alm =Ilong piece of wire anda 
thermocouple bead was made using a Dynatech Corporation 
thermocouple welder. 

¢ The other end of the thermocouple wire was connected to 
mlenadtaeadequi~sition and control unit through a junc- 


tion box. 


2. Computer Program 


MC oOcemcOonpute raptor ram |) ( FCAL) was written to 
accept the thermocouple readings through the data acquisi- 
tion system andthe bath temperature through the digital 
quartz thermometer. A listing of the computer program, TCAL 
is provided in Appendix G. This program prints all data and 
the discrepancy (i.e., the Quartz thermometer reading minus 
the thermocouple reading) as well as it stores the data on a 
computer disk. In order to convert the e.m.f. values to 
temperature, the manufacturer's conversion equation (D.1) 


was used in the computer program TCAL: 


y= ag+ a, E’ + a,E* + aE? Dea) 
meee” aa Be ae a ae E. 
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C. 


Wwiiete : 


T = temperature (°C) 
ay = 0.100860910 
a, = 25727.94369 
ay = -767345.8295 
az= 78025595.81 
ag= -9247486589 
ac= 6.97688E+11 
ae ee polo 19/2 ical 
az= 3.94078E+14 


E = thermocouple reading (volts) 


CALIBRATION PROCEDURE 


Since the temperature measurements during this investi- 
gation ranged from 25 °C to 70 °C, thermocouple fea 
bration was performed in this region. First, hot water 
(at about 75 °C) was added to the thermos flask cali- 


- bration bath. Following proper mixing, all thermocouple 


readings were recorded through the data acquisition 
system and the bath temperature was measured by the 
quartz thermometer. To obtain lower temperatures, 
small quantities of cold water were gradually added to 
the bath. 

The discrepancy (i.e., the Quartz thermometer reading 


minus the thermocouple reading) was plotted against the 


thermocouple reading as shown in Figure D.2. The 
second-order polynomial curve (D.2) shown in hve 
figure was generated using the data from alle 
thermocouples. 
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DCE 040 eon 0277043E-03 T - (ee) 
9.3259917E-05 T? 


where: 
T = thermocouple reading 
CEaoneequatton ws lay (aC ) 
DCP = discrepancy (K) 
Thus, the corrected temperature values were obtained 


Memne equation (D.3): 


m= DCP + T (D.3) 


winere: 


TL, = actual calculated temperature (°C) 
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APPENDIX E 
DATA REDUCTION PROGRAM 


[O00) STEE Wave: DRER 


1010! 
Pa 
ao 
1040 
1050 
1060 
1970 
ARO 
1090 
1100 


DATE: October iodo 
REVISED: November 1!, 1984 
BEEP 


PRINTER IS 

PRIETSUSTING “4X. “Seileeteoetion = 

PRINT USING “6X.°"O Taine data or re-processi nag prey lallcuc usm 
PRIFIT- USING “6X, °°") Plotting data = 

INPUT Ide 

[Piao HE CAlieitada 

[FP ‘bdp=1 THEN CALI SPloet 

END 

SUB Main 

COMmyeer 07> 

DIM Emf(t2),1 (12), Dlace), Daag DiatZ) , Peat2? a2) {bia eee 
DATASO. 10086091 ,257275 34963, -/5/ 34>. 2295.7 -Ueso se 
DATA -9247486589,.6.97688E+4+11,-2.66192E+13.3.94078E+14 
READ CC) 

PRs 7/01 

CLEAR 709 

BEEr 

INPUT ENTER MONTReSPAle AND TINE (MA SDD HRs) wares 
OUTED" ; Dawes 

CURBOeeta 1)" 

ENT Eels Dates 

PRINS 

PRINT 

A ehia 

Ee Month, date andatime ) Varnes 

PRINT 

PRINMSUSING “10X,°"NOIE: Progranemamensl os 

aie 

[HESENIER DISK NUMBER 2 bm 

PRINI@USING “16%, "Disk numpen (= 222 20m 

RBar 

INPUT “ENTER INPUT MODE <0=32054h2 1-Pile Sn 

IF Im=0 THEN 

peeE 

[HFUF GIVE A NAME FOR THEeRAaneDhTA-- ee 202 aires 
PRINT USING “15%. "New fie mane: 714A “D2 uf alee 
CREATE BDAT D2_file$,20 

ASSIGN @File2 TQ De file 


DUMMY FILE UNTIL Nrun KNOWN 

O1 filedS="DUMMY" 

CREATE BRAT Digi let, 28 

ASSIGN @Filel TO D1_files 

NUTFUT @Filel:Dated 

Pee - 
INFUT “GIVE AR NAME PGR THe seen Flee sea tiles 
CREATE RBDA Ie +t tlet.o 

ASSIGN slot 10 FPrepries ‘ 
Bele i 

[INFUT “ENTER NUMBER MF DEFECT Ge iGs «0-besURi alate 
LE tate 0 ier 

Ldtel=0 

Ledtes-0 

PRINT USING “16%, "No detective lGceexisce. 

END IF 
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Can a) Oy ty) a) 


meee ee ese ke CL Sena as i ae 


OO OO 2 0 0) 2) 2 WD SIN SIA ™ A NI 


870 
1880 
1890 
1900 
eae 
eZ 
cleat 
1940 


fee lidte=s| Hen 

BEC 

eT Sh BRS DErECTIVE TC EOeeatoN .bdtc | 
PRINT USING "16X,""TC 18 defective at location 
maltc. =) 

BD LF 

if Idte=2 THEN 


BEEP 

INPUT “ENTER DEFECTIVE TC LOCATIONS” .Ldtec! ,Ldtc2 
Soli eUsotNGe lon, IC are defective at locations 
END IF 

IF Idtc>2 THEN 

BEEP 

PRINTER IS 1 

BEEP 

PRINT “INVALID ENTRY" 

PEONIER IS 701 

Both 1530 

END JF 

Poe) 8F i lelt*tdtel,Ldtcz 

Im=1 option 

Bese 
Beer 
INPUT 
PRINT USING "16X%.""Uld file name: 
ASSIGN @File2 TO De files 

ENTER @Fileé:Nrun 

Bek 98 yle2:DoldS 

Beer 

PHPUT “GIVE §& NAME FOR PLOT FILE”.P files 
Seeele BUAISr filet.S 

ASSIGN @Plot IO P_filed 

molvinwoINGe 16%, IHhis data set taken on 
ENTER ®File2:Ldtel .Ldtce 

SeiwuatiGe lox.  Ihermocouples were detective at locations: ”.2(3D.4xX) : 


ae ic tes 


VE rane Mode | wed tee 


Peon Sitiowexto (ING DATA FILE .Demrr les 


rae *0e tiles 


— «140 Dells 


mame} -LdtcZ 


130 
1960 
1370 
1980 
Ue 
2000 
2010 
2020 
2030 
2040 


2120 


SS 


Cong ™! Oy Ulead 


[Sr Pa INS 


0 


i = a 


mer *ilec:1tt 

END IF 

a OUNTER IS | 

IF Im=0 THEN 

BEEP 

poltviewocllG “4%. "Select tube type" 
PRINT USING “6X.°°0=Smooth 4 inch Ref” 
SreimatG “6X. "t=Smooth 4 inch soft Cu" 
mewn “5X, 2=Smoeth 8 inch sott Cu” 
PHPUT Ltt 

IF Itt>2 THEN 

BEEP 

PRINT “INVALID ENTRY" 

GQTQ 2000 

mp) LE 

mateUT @Filel:Itt 

Breet & 

PRINTER IS 701 

ZeeteietotNG “Tax, Tube lype is: ~".0° 


Pee ent eR BUTPUT VERSION (O=(GRN. |=SHORIoS 


yet 


Di=aUrameter at thermocouple positions 
Can aie Ss. .0111125,..01143 
READ Dial) 


oa 


2200 


ae 
Ie et 


ESE 


2230 
2240 
2250 


<cb0: 
may i) 


5580 
20 
2300 


2310) 


Zac 


2330 
2340 
coo 


2360! 
2a 


2oe0 
eo 
2400 


2410! 
2420! 


2430 
2440 
2450 


2460! 
2470! 


2430 
2490 
2900 
2510 
25a 
oil 
coe, 
2% 
560 
2570 
2580 
2590 
2600 
2610 
2520 
2630 
2640 
2650 
2660 
2670 
2630 
2690 


ZA 
Zit! 


2720 
2730 
2740 
ro 
2760 
2Cee 
2780 
2730 
2600 


Dilan elt t ) 

U2=Dtameter of testmece HaSe wont tunis 
DATA .015875..01S58/sea0 

READ Nese) 

DZ=DeaCltt? 


tion to the 
15875 


Di=Inside diameter of unenhanced ends 
DAtiGaee Ol 224050127, 0127 

READ Diat*) 

P= Diet.) 


Ve=Mitside diameter oreunenvanced ends 
DATA .015875..015875,2015673 

READ Doa(*) 

Do=Doa(ltt) 


[=KeneGn ot enhanceq istic) ace 
Data. WOLG. 310 16,42022 

READ La(~*) 

Pa=wat it 


Lu=Length of tnenhanced surface at the ends 
DATA 20254, 302545 .0752 

READ Lual*) 

Lu=LuaCItt) 


Keu=Thermal Conductivity of tube 
DATA 398.344,344 

READ Keual*) 
Keu=Keuali tt) 
A=P l= Cloe2-i1 
P=PI=Do 

J= 

Sx=0 

Sy=0 

Sxs=0 

Sxy=0 


2) 7/4 


Repeat: ! 


IF Im=0 THEN 
HVE e eo RECOVER 2o30 


PRUNES | 

PRINT USING 7 SELECT GHiet lt) tam 

PRINT USING he 0 OIE Ue eee 
‘T=SE] ae FLUX 


Pe Wit 99K 
ESCAPE” 


rE Ke 


PRINT USING ON =e eed = 
PRINT USING “4x, “NOTE: Re ‘9 = 
BEER 

INE UES 

IF Ido=0 THEN 3680 


LOGR? Wiese ent Flux 
IF Ido=1 THEN 
BUTPUT SZ Ose AR AP Se AGS SIRS 


BEEP 

INPUTS ETERS DESIRED: Qdp Dade 
PRINY USING 4x" DEST REDS0ap 
Eres 1000 

POR ae tee 

OUT EU a Steel 


Sum=(0 


ACTUAL Qdp 


a 


POR Ji=! 10 5 

BNIER 7/0976 

Sum=Sumt& 

WEA T Jt 

IF I=! THEN Vol t=Sum*5 

IF [T=2 THEN Amp=£— 

NEXT I 
Aqdp=Volt*Amp/(PI*=D2«L) 

IF ABS(Aqdp-Dadp)>Err THEN 
IF Aqdp>Dadp THEN 

BEEP 4000,.2 

BEEP 40900, 
BEEP 4000, 
ece 

beer 250..2 
peer 2ou..2 

Por col..e 

END IF 

PRINT USING "4X ,MZ.3DE.2X.NZ.3DE":Dadpe.fAaqdp 
WAIT 2 

GOTO 2780 

Bese 

peer 

ee USING 4xX9MZ,.3DE .2X.4Z.SDE° :Dadp.Agqdp 
Err=900 

WAIT 2 

BuT0 2780 

aL 

END IF 


nN 
Ne 
r 
ee 


' LOGP TO SET Tsat 


IF Ido=2 THEN 

AUTPUT 709;:"AR AF33 AL33 VRS" 
ONUTPUT 709:"AS SA" 

BEEP 

INPUT “ENTER DESIRED Tsat”’.Dtld 
PRINT USING “4X.""D Isat Amisae 
Old=0 

alee 7 UIE lta 
Atld=FNIvsv(Eliq) 
Beae=PNPsat(Atld) 

IF ABS(Atld-Dtld)>.2 
IF Atld>Dtld THEN 
feer 4000,.2 

BeeP 4000,.2 

peer 4000..2 


BEOL 
pee 250 ..2 
olen 


BEEP 

Meee 250... 2 

eA LE 

err-Atid-Uld 

Old=Atld 

Sern USING “4X ,.4¢C4DD.00, 3X7" :Dtld.Atid.err. Psat 
WAIT 2 

SOTO 31390 

Be 
tewAbotntliGg-Dtld)>. | 
IF Atld>Dtld THEN 
BEEP 3000,.2 

Beer oU0U. 2 


Diff ea 


THEN 


- 


THEN 


oS 


3420 
3430 
3440 
3450 
3460 
32470 
3480 
3490 
3500 
3510 
3520 
3520 
3540 
3 0) 
3560 
3570 
3580 
3590 
3600 
3610! 
3620 
3630 
3640 
3650 


3660! 
56a0 ! 


3680 
3590 
3700 
3710 
3720 
3730 
3740 
3750 
3760 
3 a 
3780 
3/790 
3800 
3810 
3820 
3830 
3840 
3850 
3860 
3870 
2880 
3890 
2900 
33g 
Soe 
3330 
3940! 
Bais.) 
3960 
aoa 
3989 
3990 
4000 
4010 
4020 


ELS 

REEF VEO. 2 

BEEPS Wie. 

END Le 

Etre -ulld 

Old=Atld 

He USING...“ 4X24 NDDRD Ss eet lid At lane nee oot 
ALhee 

GOTO 3190 

Piece 

Geer 

ETre-felaan ld 

OQOld=Atid 

PRINT USING "4X  4¢(4DRSpUes “Dtlod.Gtiawerr-r san 
WAIT 2 

GOTO" 3190 

END IF 

END IF 

END. IF 

ERROR TRAP FOR Ido QUT OF BOUNDS 
IF Ido>2 THEN 

BEEP 

GOTQ 2620 

END IF 


TAKE DATA IF Im=0 LOOP 

BEEF 

INPUT “ENTER BULK DILSG store 
OUTPUT 709:"AR AF25 AL36 VRS" 
FOR T=t 10 t2 

QUT EUs: AS SA” 

Sium=0 
PORSGh = eel 20 
Ee ee 


Sum=Sumte 


NEXT Ji 

Emf (1) =Sum/20 

NEXT I 

OUTPUT 709;"AR AFG? ALG3 YRS" 
Fieel—ies bl 2 

BUTPUT 709;:"AS SA" 

Sum=0 

EAR Jt 10 20 

ENTER 709:E 

Sum=Sumte 

NEXT Ji 

IF I=! THEN Vr=Sum/20 

IF [T=2 THEN Ir=Sum/20 

MEXT I 

ELSE 

ENTERS oF lee Bop, folds ,fmt ( aeey ral 
END IF 

CONVERT emf’°S TO TEMP ,YOLT,CURRENT 
Twa=0 

EGR Ol = ap Bee 


IF Idtc>N THEN 

LF I=lditeiv GR l=Ldtec THEN 
iho == oso 

GOTO 4060 

END IF 


94 


4290! 
4300! 


4310 

320 
4330 
4340 
4359 
42360 
4370 
4380 
4330 
4400 


4410! 
4420! 


4430 
4440 
4450 
4460 


4470! 
4430! 


4490 
4500 
4510 
4520 
4530 
4540 
4550 
4560 
4570 
4580 
4590 
4400 


4610! 
4620! 


4630 


eNy IF 

TCI =FNTusulEmf (1)? 
Pee <S IREN itia=IwasT< 1) 
NEXT I 

itva= Tuas (S-=ldti? 
igeai= 1 (2) 
Me=C(C10)4IC 1172) 
Tsump=7T¢(12) 

Amp=Ir 

VYolt=Vr<25 
Q=Volt«Amp 

IF Itt=0 THEN 
Kou=FNKcu’ Twa) 

Boe 

meu-Keual ltt? 

END IF 


FOURIER CONDUCTION EQUATION WITH CONTACT RESISTANCE NEGLECTED 
Tw=Twa-Q=LOG(D2/D1)/¢(2=PI*Kcurl? 

Thetab=Tw-Tid 

IF Thetab<O THEN 

BEEP 
INPUT ““TWALL<TSAT (O=CONTINUE. 
IF Tev=0 THEN GOTO 2590 

IF Iev=! THEN 5030 

END IF 


COMPUTE VARIOUS PROPERTIES 
le tim=eNifilm¢ lw, tld) 
Rho=FNRho( Tf ilm) 

Mu=FNMut( T£ilm) 
K=FNK(Tf£ilm) 
Co=FNCp(TFilm) 
Beta=FNBeta(Tfilm) 
Ni=Mu/Rho 

Alpha=K/(Rho*Cp) : 
Pr=Ni/Alpha 
Psat=FNPsat(Tld) 


t=EHD)” . lev 


COMPUTE HEAT TRANSFER COEFFICIENT 

Hbar=190 

me hbar@~F/(Kcou*A))~ .5*Lu 

Tanh=FNTanh(Fe) 

Theta=/lhetab«fanh/Fe 

PRINTER IS 701 

Pee lhG “4%. 7CiX ~MZ.SUE? :Hbar.Fe, lanh.Thetab.Theta.Beta,Ni 
Kx=(9.81=Beta*JThetab*«Do 2=Tanh/(Fe*Ni*Alpha)) .165667 

Pelt (soo3/Fr) ~~ (9/16)) «8/27> 
foare=;/Do~( .6+.38/*xXx/Vy) 2 

ieenes “ibar-Hbharc)/Hbarc?)>.00! 
fear -(t4barthbarc)*.5 

GOTO 4440 

END IF 
Ql=(Hbar*P*Kcu*A?) 
Pe-W—2 ~G)| 
As=PI+«D2*L 
Qdp=Qc/As 
Htube=Qdp/Thetab 


RECORD TIME OF DATA TAKING 
IF Im=0 THEN 


THEN 


Fone tabs latin 


dis 


4640 
4650 
4660 


4670! 
4630! 


4690 
4/00 
4710 
47290 


AQT sae 


4730 
4740 


4750 


PRINT 


Oe 0S 


Pelee 


Emcee |olas 


ED a 


DUTPUT DATA 
PRUNE R SIS 
Pech 


To PRONE 
701 


IF lTov=0 THEN 


PRINT USING 
Bop. Told$ 
PRINT 

USING 


PRINT USING 


().. 038 


4760 
4770 


etab, Ivyarsat.lsump 


4780 
47930 
4800 
4310 
4820 
483 

4340 
4850 
etab 
AGB 
4870 
4380 
4890 
4900 
43140 
4920 
4930 
4940 
44950 
4960 
a37 i 
AI80 
4990 
5000 
50g 
Suu 
S30 
5040 
Bat 


5060 


Pr LPT 
PRINT 


USING 
USING 


PRINT USING 
PRINT USING 
END IF 
IF Iov=] 
IF J=! THEN 
PRINT USING 
END are 
PRINT USING 


END IF 


"t1X," “Data Set Number = "",DDD.2X%.”°" Bulk Wil % =)" pee ee 


TX ee GaNos | Z 3 a 9 6 / 


“EEX. Temp 3°". BCtX MDD.DDI":1T¢1),. 762). 103) . 1 eee 


IT el Maed Thetab [Caer Psat Tstines 
"11X,2(HDD.DD.1X) .MZ.DDDE,.1X,2¢41X,.MDD.DD),.2%.MDD.D  : lite eee 


a9 sooges 


X 
1tX,°" Htube Odp 
“11% ,4MZ.3DE, xtlz SUE” (tube. Gar 


THEN 


"Lika SRUN SNe ls ack Htube Qdp Thetab = 


"t2X .DDD.4X,DD.2X.MDD.DD,3¢1X,MZ.DDDE)”’ :J.8o0p. I1d, Httise, Ga 


on OMT 


Beet 


Ue eer STORE ES sDe in 


END 
ea) ema ate 
ee eer 


IF Im=1 GR Ok=1 


SET ¢(t=Y.0=N) ?".Qk 
Im=! THEN J=J+) 

Im=0 THEN QUTPUT GFilel:Boo. Told3,.Emf (=) ,Vr.Ir 
THEN QUTPUT @®Plot:Qdp, Thetahb 


IF Im=0 THEN . 


Bee 

PUR wie 
Nrun=J 

[FE “Gewen<> | 
IF Go_on=1 
ELSE 

IF J<Nrunt] 
BD ie 


THERE BE ANBTHER RUN ()—) 30S 7. Gowen 


THEN S030 
THEN Repeat 


THEN Repeat 


Me le lS, 


BEBE 
Peri 
PRINT USING 


1. D2 Ai ele 


aU70 
9980 


ASSIGN 


File 


“TGX. "NOTE: °".2Z,.°° data rims were stored int 1)e ee re 


al 


Ce ee Vina | 


ASSIGN @Filel 


ENTER @File 


TO Dt files 
Wate oy tic 15 eee ele te 


GUTPUTB@Rile2 Dates jadiget lateo. bee 


FOUrelS4 
ENTER File 


Ne oie 


Leb opal aldS> Emi <oes\ nlr 


OUTRUT GE ilez- Bop. loldS.emt (ovr ir 


NEXT I[ 


ASSIGN @Filel 


0 


PURGE. UR = 


ENO Se 


a6 


al BEET 
SZUU PRIN | 


Seamer RINT USING “10X.° “NOTE: °°" .ZZ, °° X=-Y pairs were stored in plot data file 


Seeerufn :J-1,.P files 

eeu ASSIGN @File2e 10 ~ 

wee fAsSIGN @Piat TO 

eee 6M BEEP 

5250 INPUT “LIKE TO PLOT DATA (1=Y,0=N)?",Ok 
5260 IF Ok=1 THEN 

me =CALL Plot 

eee) END IF 

5290 SUBEND 

5300! 

Soo! CURVE FITS OF PROPERTY FUNCTIONS 
me) DEF FNKcuCT) 

pee’ OFHC COPPER 250 TO 300 K 

S340 Th=T+273.15 Vea K 

pe «=k = 434-,112*Tk 

S60 RETURN K 

pee FNEND 

5380 DEF FNMucT) 

eeee') 170 TO 3260 K CURVE FIT OF VISCOUSITY 
page Tk=1+2/75.15 ‘Cc 10 K 

S410 Nu=EXP(-4.4636+¢(1011.47/Tk))>~1.0E-3 
S420 RETURN Mu. 

Sas FNEND 

[440 DEF FNCp(T) 

Sou is0 10 400 K CURVE FIT OF Cp 

5460 Tk=T+273.15 iC are 


eeepc e=. 4188+! .65007E-3*Tk+1.S1494E-6*Tk 2-6.6/853E-1U0=Tk 3 


9480 RETURN Cp 

5490 FREND 

See VEF FNRhoCT? 
weapoeesik=1+2/3.15 nee TK 

meee A-1-(1.8*[k/753.95) IK TOR 


Soo) Ro=36.32+61.146414*X°(1/3)+16.418015<X+17.4/76838=XxX 


5540 RETURN Ro 

eee) EMEND 

5560 DEF FNPr(T) 

5570 Pr=FNCp¢T) =«FNMuCT)/FNKCT) 
S580 RETURN Pr 

5590 ‘FFNEND 

See DEF FNKCT) 

5510! T<360 K WITH T INC 
S620 K=.071-.000261+T 
Be30 RETURN K 

modu ) FNEND 

5650 DEF FNTanh(Xx) 
S660 P=EXPC(X) 

moe Q=EXPC(-X) 

S680 Tanh=(P-Q)/(P+Q) 
56390 RETURN Tanh 

5700 FNEND 

mee, DEF FNivsvcV) 
meee COM /Cc/ CC7) 
ey7ou 6 =C(0) 

S740 FOR I=1 T0 7 
57590 T=T+CC])«V I 
5760 NEXT I 


em Oo. 7 4e7°9S4E-24+1*(9,.027/043E-3-1*(-9.3259917E-5)) 


a7 


eB ZO Xk —e 


5780 
a7 
5800 
on 
5820 
5830 
5840 
5850 
53860 
5870 
5880 
S570 
5300 


ote 


S920 
3430 
5240 
S3a0 
Soe 
aay 0 
5380 
S250 
6000 
6010 
6N20 
6030 
5040 
6050 
5060 
6070 
£080 
6090 
5100 
6110 
5120 
6130 
5140 
6150 
5160 
6170 
65130 
6190 
5200 
6210 

220 
5280 
6240 
BZo0 
6260 
Ber 
6280 
b290 
6300 
5310 
6320 
5330 
6340 
53a5n 
6360 
6370 
6380 


RETURN. 7 

FNEND 

DEreritbeta [> 

ROO hemow ht. |?) 

Rom=FNRho(T-.1) 

Beta=-2/(RoptRom)=*(Rop-Rom)/.2 

RETURN Beta 

FNEND 

DEP ee Ni eelm lu, bidd 

TFEilm=€Twt+Tld)/2 

RETURN Tfilm 

FNEND 

DEF FNPsat(Tc) 

0 T0230 deg F CURVE Bi Supe scat 

Tt 26=1e+32 

Pa=5 .9455925+TF=*¢€. 153952082+TF#*(1,. 4940963E-3+7TF*9.6G150671E-6)) 
Poe Fae e./ 

IF Pq>O THEN ! 
Psat=Pq 

LAs 
Peagt-rpaecs.,9-/!4.7 
ENP Ie 

RE RURINer oat 

FNEND 

SUB ME lot 

DIM SeG3p 

PTE Gee ly 

PRINTER SiS) 


+ SG arg 


Idv=1 
Beer 
INPUT “LIKE DEFAULT VALUES FOR PLOT (1=Y,.Q=N)?",Idv 
Bee 
PRINTUUSING “4X. “Select Mption: 
PRN ieUS ING)" 4X 2 Gee ec teeeee lee ies 
PRINT USING “4X. (he versus dees 
PR eSNG OX oe Ze ee 
PRINT USING “4X,°"3S “h=Tatio verstis delta) 
INPUT Opo 
IF Opo=3 THEN 
Cea 
INPUT “SELECT TUBE DIAMETER (0=.75,.1=!.0 IN)” ,.Itd 
END IF 
PR eels 7 (lo 
IF Idv<é>1 THEN 
Bee 
ee “ENTER NUMBER QF CYCLES FOR X-AXIS".Cx 
ceo “ENTER NUMBER OF CYCLES FOR Y-AXIS" .Cy 
EP 
INPUT “EMFER MIN X-VALUE “MULTIPLE Gee a> niin 
BEEP 
LVPUTS CENTER WIN y—VALUE “(MUIR PrE he I meen 
ELSE 
Cy=3 
IF Opo=0 THEN 
Cx=2 
Xmin=1 
Ymin=!1000 
END IF 
IF Opo=!t THEN 
Cx=2 


ae 


SoU 
6400 
6410 
6420 
6430 
6440 
6450 
6450 
£470 
6480 
5490 
6500 
sya 11) 
6520 
6530 
6540 
6550 
6560 
6570 
6580 
boJ0 
6600 
bE10 
6620 
5630 
6640 
£650 
6660 
6670 
6680 
6690 
6700 
5710 
6720 
6730 
5740 
B67 aU 
6760 
6770 
6780 
6790 
6800 
6810 
6820 
6830 
6840 
6850 
5860 
6270 
5880 
6830 
5900 
Erne 
6920 
6930 
6940 
6350 
5460 
63970 
6330 


Amin=| 
Ymin=100 

END IF 

IF Opo=2 THEN 
Cx= 


Xmin=1000 
Ymnin=100 

END IF 

tte tipe=3 THEN 
Cxr=2 
Cy=3 
Xmin=] 
Ymin= | 
END IF 
END IF 
Beer 
ee Leh 
moet Se OF 100,0.100;TL 2.0;:" 
Sfx=100/Cx 

eey-100/Cy 

aoe Ew 0.0 PD 

Nn=9 
FOR I=! TO Cx+l 
nace xXmine 10) €1—))) 
IF J=Cx+!t THEN Nn=!I 
FOR J=1 TO Nn 

IF J=! THEN PRINT 
IF J=2 THEN PRINT 
Xa=Xat*J 
X=LGT(Xa/Xmin)*Sf x 
ooo PAX, 0: 
NEXT J - 

NEXT I 

fer FA POO, sPU:” 
coom)6 6 6PhU UPA O.O UPD” 


Nn=93 
Giles 


ieee 
ce 


ee 


FOR [=] 
mwat—imin= i (1-1) 

IF [I=Cy+! THEN Nn=l 

awR J=1| TO Nn 

Sai pHen PRINT “TL 2 0° 
eee -2 THEN PRINT “TL 1 0” 
Pa-Tat~J 
Y=_LGT(Ya/Y¥min) *Sfy 

aa FAO, SY. YT" 


NEXT J 

NEXT I 

eet. FA 100 IL 0 2° 
Nn=39 

aie l=| Th Cxti 
Rat-Amin*10 ([-1) 

IF [T=Cx+! THEN Nn=l 

FOR J=1 TO Nn 

Peete thEN PRINT “TL O 2 
IF J>1 THEN PRINT “TL O 1” 
Ra=Xat=J 


X=LGT(Xa/Xmin)*Sf x 
ieee SX, | 00: 
Hex) 2 
NEXT I 


XT" 


99 


inom mor 2500, 1800, 8300.6900; — 


90) 
7000 
7010 
7020 
7030 
7040 
AUS 
7060 
7070 
7080 
7990 
7100 
7110 
7120 
7130 
7140 
7 lou 
7160 
770 
7180 
70 
A200 
7 cae 
7220 
7ae0 
7240 
7290 
7260 
7270 
7280 
Caan 
(eu 
7 a0) 
7320 
(350 
7340 
7350 
7360 
7370 
7380 
( Saw 
7400 
7410 
7420 
7430 
7440 
7450 
7460 
17 0 
7480 
7490 
75-0 
731.0 
7s20) 
780 
7540 
( so) 
7560 
7570 
7360 


PRINT “FA 100.100) PUDPaR oie ier pr 
Nn=9 

FOR I=! TO Cytt 

Tac =e l—1) 

IF [T=Cy+? THEN Nne=l 

FOR J=1 TO Nn 

IF J=! THEN PRINT "TL O 2" 

LE JO) THENSPRINI “Tle Ge 
Ya=Yat*J 


Ve EG ora Vmin) at y 


ER ENG PDA I00 Ye 
NEXT J 

NEXT I 

PRINT “PA 100,100 PU" 
PRINT “PA 0,-2 SR 1.5.2” 
Tishaa ne wm) 

FOR IT=1 TO Cx+t 


hae sien ee 1 
X= LGT(Xa/Xmin) *Stx 
Seer eee et A yore 
IF Ti>=0 THEN PRINT 
eC LN 


"CP =2.72 (EERO RR =2o72 Picea een 
"CP 24-27 LBd0- PR U2 IED ele eee 


lor= ie t 

NEXT I 

Polhieee wee rn O,0* 
T1=LGT¢CYmin) 
Y1Q=19 

POR [=a 10) Cy+i 


a= (neo 1 — | 
Y=LGTCYa/Y¥min)*Sfy 


PR SPA ey 

PRINT "CP i. = 2 eB Gee eee ge le al tere aes 
Ti=li+1 : 
Ne cle 

IF Idv<>t THEN 

BEEP 

INPUT “ENTER X-LABEL".XlabelS 

eee 

INPUT "ENTER Y-LABEL”.Ylabel$ 

ase 

IF Opo=0 THEN 

Piabels= Tw-lsat Ge 

Voabemo= auth /m 2)" 

END IF 

[IF Opo=! THEN 

AlabeolS=" {n= leat. ch 

YlabelS= bh CA/mecek 

END IF 

IF Opo=2 THEN 

lai | sa= cacy mae. 

TlabelS= hech/m c.k)) 

END IF 

IF Opo=3 THEN 

Alabelh="Tu-lsat (Fh) 
YlabelF="hCenh)/h(smooth?)" 

BNUe he 

tty ae 

Pe Rl Sxte we i ePS 20. -{6 CP’ :-=LENCX label a)7 2aaemee 
PRINT “PA =14.50 CP 0.°:-LEN CY label$)7 2-576. Dees 


PRINT? "CP Ge) 1 


75390 Repeat:! 


100 


‘XlLabel$: 
lane) hae 


eee 


fbUU 
7610 
7620 
7630 
7540 
7650 
7660 
oe) 
7680 
7630 
7700 
7710 
7720 
7730 
7740 
7750 
7760 
at 
7780 
Ba) 
7300 
7810 
7$20 


WNDU & WIT — 
ey eee Gee ee ey 


BEE 


INPUT “WANT TO PLOT DATA FROM A FILE (1=Y,0=N)?" Qk 
IF Ok=1 THEN 
Beer 
INPUT “ENTER THE NAME DOF THE DATA FILE’ .D_files 
ASSIGN @File TQ D_filesd 
PEEP 
beer 
Bee "ENTER THE BEGINNING RUN NUMBER" ,Md 
ae 
INPUT "ENTER THE NUMBER OF X-Y PAIRS STORED" .Npairs 
BEEP 
PRINTER IS | 
Berit USING 4x. "Select a symbol?” 
PRINT USING Pie | (star Sam ees. Sain 
PRINT USING "aX, hn oeecire le 4 ee ras | 
PRINT USING "4X%.""S Rombus”"” 
Bre iSiNG 4X2" 6 Right-Sside-up triangle” 
PRotimlictINnG 445 7 Wse-side-down triangle” “ 
INPUT Sym 
PRINTER IS 705 
Bein eu DI 
IF Sym=1 THEN PRIMI “SM*" 
Ce sie 2 ele eect Sven 
[eeeem=3 (HEN BRINT “SMo” 
IF Md>{t THEN 
FOR J=1 TO (Md-1) 
ENTER @File:Ya.Xa 
NEXT I 
END IF 
FUR t=! TH Noairs 
aeR OPwhe:ya,Xa 
IF Qpo=1 THEN Ya=Ya/Xa 
IF Dpo=2 THEN 
W=Ta 
ha-ida/ Xa 
Xa=f) 
END IF 
IF Bpo=3 THEN Ya=Ya/FNHsmooth( Xa, [td) 
X=LGT(Xa/Xmin) =Sfx 
Y=LGT(CYa/Ymin) *Sfy 
femovmes THEN PRINT “SM” 
IF Sym<4 THEN PRINT “SR 1.4,2.4" 
Geo “PA .X.Y.°" 
feoymes THEN PRINT “SR 1. ae 
IF Sym=4 THEN PRINT "“UC2. és. 7 eae r a eee Q: 
feesyosoelheh PRINT “UCS,0.39.-3.-5,-3.6. eee 
ieeovymeo THEN PRINT “UCO.5.32,939. 3, -8,-6, 0. . " 
IF Sym=7 THEN PRINT eI ee oie a 
NEXT I 
Beer 
ASSIGN @File TQ * 
GOTO 7600 
aie Ee 
PRINT “PU SM” 
Pee 
INPUT “WANT TO PLOT A POLYNOMIAL ¢€1=Y,0=N)?" ,Go_on 
IF Go_on=1 THEN 
Pee 
INPUT “ENTER LOWER AND UPPER X-LIMITS".X11,.Xlu 


LO 


8200 
9210 
8220 
Beg) 
8240 
8250 
8260 
8270 
8280 
8290 
8300 
Sow 
8320 
3330 
3340 
8350 
8360 
8370 
3380 
e290 
8400 
8410 
3420 
8430 
$440 
8450 
3460 
8470 
93480 
8490 
8500 
8510 
B20 
8530 
$540 
So a0) 
e560 


FOR Xx=0 TO Cx STEP Cx/200 
Xa=Xminx!0 Xx 

[F XaSxXl? OR XaoxX liu THE eeu 
Ya=FNPoly (Xa) 
Y=LGT(Ya/Ymin) «Sty 
X=LGT(Xa/Xmin)*Sfx 

IF Y<0O THEN Y=0 

le yoo THEN GOT s220 

PR TN ie en oe 

NEXT Xx 

END IF 

PRINT “PU PA 0.0 SPQ" 
SUBEND 

DEF FANHSmooth(X, Ltd) 
Hs=FNPoly(X)/X 

IF Itd=1 THEN Hs=Hs*.83347 
RETURN Hs 

FNEND 

DEF FNPoly (xX) 


Poly=-4) 4123713E#2-X=(6.. 31239178 +2-%*3 7 Wiech tae 


RETURN Poly 

FNEND 

DEF FNPyvst(ITsteam) 
DIM K¢8) 


DATA ~-/.691234564,-26.08023696.-168.1706546 .64, 23285504 , “lie. 9e4Gae 


DATARS I? 32,20. 97 So Gy on lee ain 


READ K(*) 

l= Cleteamte 7c. 15/64/7435 
Sum=G 

FOR N=0 TO 4 
Sum=Sum+K (CN) *C1-T) 7 ON+4+1) 
NEXT N 


Br=Sum/(CT=C1+K(5)*(1-] 04K CG) 41-1) 2)9-(€l—-1)7 (hee oe 


Pr=EXe Car» 
P=22120000*Pr 
RE TUR Ieee 
FNEND 


Tbe 


APPENDIX F 
AN EXAMPLE OF REPRESENTATIVE DATA RUN 


fonthemeacerwand time “2:0 t216 25:24 


NOTE See todram Tanne DRP 
Dick onttmber = HP 
Ald file name: WHOS 
Livtieeedeea Set “taken on pie ii 45729 
iferm@courles were Getecriwwe at) locations: (I 
Tube Type 13 EB 
Data Set Number = elk eee. = 051) ere el eae, 
Te No : 1 ae 3 A : F, / ¢ 
Temp Beiee So.ly thh.04 7f?3 78.58 6.22 76.21 59.24 
Tua ielerae Thetab (Vann Se oa Toump 
6.34 10.03 %S.S6NE+01 AW aah 1380 
Htube Qdp 
r.S45E*03 8.5886 +04 
Nata Set Number = ecules =. OG flee | tas 52 
IG Noe 1 é 3 4 S 5 i Fr 
Tempo eee eee oo) Ue doe oem omen. co tea 97.76) | 647 04 
Twa Mist ad Phe tab Sete Sant Teump 
49,84 19.07 3.939E+01 10.230 Ren = 125,49 
Htube Qdp 
T.ES2E+03 8 4.538E+04 
Data Set Number = Swerik Oeies. = 020 Deets te 3G. (7 
re Nee 1 a 3 4 a f ? €? 
Tenp Co ees oe eee ce ey eel Soemon. 0 932.96 
Vwa_ livcad The tat ivapie orsat Tsump 
foe 00) 6. SO 7E+01 ee amet =| beeS 
Htube Wdp 
Pe4Z0E+Ge 1,7496+04 
Data Set Humber = Aoi Uae. = 1g eee 1 249 Sts 
TC Na: 1 2 3 4 i f- 7 g 
Temp See ie re re ese dees oe, 79.00 2b. 
Twa ine hetae [Vania ee sat Tsump 
68.79 10.07 1.S5536+01 ele 3.988 -11.9 
Htube Qdp 
peewee o6..o64E 403 
Data Set Number = Reet pomlet lne me wane bain Mieethes thas? OZ 
reno: ] a 3 | S; 6 7 $ 
Temp ee et eG terte Selatan ss § PFT! 1.659 
rare ID isfa [he Gab a ae eee Sle Tsump 
c1.988 9.95 1!,1896+0T eed? Reet wei. | 
Htube dp 
foe te 863. GA0E tis 
Nata Set Number = Be sulk, Unie. = 90 0 Pieter te 5246 
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me ny A 


TC Neo: ] = zh z 
Temp : (8.14 7,58 17,82 ti7es7? 19.18 5812.03 8.7 
Twa i lsinenal Thetab |Waert sat Tstump 

17.35 ~ 10.19 7 eee ioe Ot) 3396. Jae 

Htube Qdp 

| 72oe Oe” “esse. 
Data Set Number = 7 Veale Oi) 7 ee luiwali2el2 | ae 
[ee 1 ec 3 4 5 6 7 e 
Teme 15,58 {5.08 15.14 1t4,.87 15.86 15.04 1S 4p ee Se 

Twa Tliad tie ba Tvapr a Tsump 

S22 10 02 Saal Lewes 3.390 Vee 

Htube Qdp 

IP22(B+02 fae ae 


NOTE: O7 X-Y pairs were stored itn plot data file POS 
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110 
110 
Pou 
129 
141) 
of 
1690 
170 
130 
190 
21h!) 
a0 
2c) 
aoU 
2 4f) 
24] 
242 
244 
245 


¢ 


APPENDIX G 
LISTING OF CALIBRATION COMPUTER PROGRAM (TCAL) 


Belle SHANG aires 


meee foe). October 4, 1984 


Gall /Ce/ C7) 
Dati emt (475,18 4s) C4) 
Commodi af c/n gobo. -6m7 sone 39,7 002001060 
DATA -3247486589.6.97688E11,-2.66192E13.3.940/78E14 
READ CC(*) 
CLEAR 709 
ce ales 
aT “ENTER MONTH. DATE AND TINE (MN: D0:HH:MN:SS)”,8S 
OUTPUT 709:"TD":BS 
COE ioe 7 De 
me 7 0 erie 
ar 
Cee NAN VASHAR DE CUPY Clay, 0=tl) 2? ,I1hp 
IF Jhp=1 THEN PRINTER IS 701 
[If Ihp=0 THEN PRINTER [5 |] 
eS UStiGawite Tantaemgate and time = "",14A 1AS 
(ieGT. ENTER @INPUl HOVE <l=s054A, 2=FILE)™.Im 
IF Im=1 THEN 
Bile 
INPUT “GIVE A NAME FOR DATA FILE",D_files 
CeenlE BOA DL files, 20 
clase 
Beles 
pe SOV Ewe lh SxS iiG Pile yOat ites 
INPUT “ENTER NUMBER ODF DATA RUNS STORED" .Nrun 
REE 


EEP 
INPUT “GIVE A NAME FOR OUTPUT FILE" Of iles 
CREATE BDAT OF i led, 19 
ASSIGN @ ileo 10 Ofiles 
K=] 
END IF 
ASSIGN @ ile TO D_files 
IF Im=] THEN 
Bt.EP 
iets SNIER BATH TERPERATURE {lobath 
Beebevtac ios TeRZE 
WAIT 2 
Bene? 133 io bath 
OUTPUT 709;:"AR AF29 AL33 VRI” 
FOR [=0 [QO 4 
eur 704: AS SA” 
Pelee, Ootvemt (1) 
WeExT I 
alee 
Biteneslile?| bath, emt t=} 
END IF 
J=Jt+] 
pone =0 olf 4 
TCI) =FNTvwevi ABS (Emf ¢(1))) 
Wel | bath-T¢]) 


NEXT I 
= eet 
FOR L=0 TO 4 


OUTPUT @Fileo: TCL).D0C) 
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S45 
548 
549 
Sat) 
S60 


620 
£30 


eyo mae 
GOTO 700 
ETD ali 
PRINT 


PRINT SUSING “Wx, Ba theevemperature i ec) (2) a 
PRINT USING “10X."" Thermocouple readings (Deg C):""" 


PRINT WSINGs 10%, 20302 DDR a ieee sx sc 

PRINT USING “l0X,°"Diserepancies (Deg Gp | 
PRINT USTHG. “LEX. SCMZ  Ditteixe 245 a) Ge 

IF [m=Z2a0he 700 

BEE 

INPUT “OR [0 ACCEPT © THIS SE! l= iG-N) 7 sae 
IF Oks=0 THEN 

Jedi 

GATN 390 

bls 

OUTPUT @rle;sT_bath,cmf (+) 

END JF 

Beer 

he ime tee 

INPUT “HILL THERE BE ANOTHER SET < l=) (Sth = Gamo 
IF Go_on=1 THEN 399 

ailsls 

LF. ale non Heh 330 

END I[F 

FR LT 

IF Im=1 THEN 


PRINT USING “l0X.""“NOIES  “.DG.  dalvaese ts -aremetranea 


Eieoe 

PRINT WSING “10%, ° "NOTE: Above analysis was Cerierpmed | + one 
PRIPT 40SING 1k 2) = Output data are stared in file 

ties 

SCLC!) Fale [Oi 

END 


OFF FNIwsvcEmf ) 
COR Wey wee?) 
Sum=C(0) 

FAR LI=1 TO ? 
Sum=SumtCCL)<«Emf | 
NEXT | 

RETURN Stam 

FMEND 
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(Deg Cy) ° ; [aber 


in file “1 4A 


iece : 14A":D_ 


"(4A Opi Tes 


AR Eee 
SAMPLE CALCULATION 


Data run number 5 (Saturation temperature was 10 °C and 
heat flux was about 20 kW/m*) was chosen for the sample 


eeeeculation. 


PeeeelGol-SECTION DIMENSIONS 


Ds = 0.01588 (m) 
15: = 0, 0270 (in) 
Di = 0.01143 (m) 
DS = 0.01588 (m) 
L = 0.20320 (m) 
i, = 0.07620 (m) 


Pee SASURED PARAMETERS 


ye = 350) CVO les) 
ie = 221 ACV lela s@) 
ae = a ae) 
> = 305 25 (-C 
T 3 eee Oi ee) 
T 4 == 33449 (°C) 
TS eee ee) 
ihc = 31.38 (°C) 
> =e) CoO Ole) 
T 9g = Seno © () 
ee = EO: 00 (2G) 
ee = See Omen tke 
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OUTER WALL TEMPERATURE OF THE BOTeE iG TUS 


=T D = 3.1416 ino 530 6 Ome 4am) 

= 7 (D2,- D2.)/4 = 351416 x (OvOmoas* - 00 
= 9714 x. 10 = iGne) 

=" 25V. = 25 KNOG = 89 0 Gare ) 

= ie 2. cain) 

= “VI =969 5O0MEe2 27 = 20a Gre 


Oo Hes ee UC 
| 


8 
lav¢é ( i) pogor =P 8 «/ 2. CAG 


n=] 
Tot Tyg Qy( INn(@,/ D,)/2 TL k,) 


To =)06|633.72 - 203.17(1n(15.88/11.43)7 2x3. 1416x203 eee 


eee 7 Ce) 
Ty - iT 6 
a 2 = _ = Oo 
6 = Tp- 1.5 33.57 = 10q00s 23a) 


PROPERTIES OF R-114 AT FILM TEMPERATURE {REF. 24} Rie 
2 


Tas Cee he Oe = 21.78 (ec) 

log U(10 *Ns/m*) = -4.4636 + 1011.47/T; 

Lb = exp(-4.4636 + (1011.47/2345 73) oeoeoe 

Ve Se) oe TO ON itn = 

P = 36.32 + 61.14 3: °°? + 1GUze eee 
ee CC lpg, tem 

J. °=) - /T, 

T, = Gritical Temperature) see (R) 

T . = 294978 (K) = S53 Choe 

j = 1] 530760) 753295) - 3022 
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= Seo sommelby £> eo 49 (ke/m ) 

ees Ey Pre 2 428 10.7 m/s) 

nw /mok = O07 10 =O 000ZE), Te (Tp in °C) 

ke = 0.0653 (W/m.K) 

C (kJ/kg .K) Or come O Let elo le x lO Te = 

poooex 10 a eran K ) 

G, =) ECOZG kek) 

eee ty ee Ga ee x 105 Cre) Ss) 

mee = -p/ &@7)/ 2 

B = -(1465.95 - 1466.49) / ((294.98 - 294.78) x 1466) 
Bee =- 0.00184 (1/K) 
meee = V/® = 5.464 


Pere AT-FLUX CALCULATION 


Average natural-convection heat-transfer coefficient 


meenen-boiling ends: 








— ce 
l 
gi Dse@,Tanh| cs ) reel a 
Oo a ; 
HP ,1/2 
: eh ag 
h = DT 0.60 +0.38/7 D716, 872 


O i Or oye) 


hee? (Wim? .K) 


Heat-transfer rate through non-boiling ends: 


Or 


= 7 / 
Q = (hp k,A.)*’? 8 Tanh(mL, ) 
m = (h p / kK-A y=) 23504 ee mn 


Q. 13.15 (W) 
0 es 2% 00 = 7a oreo (W) 


Heat flux through active boiling surface: 


Q = Oe OE ads 203.17 = 26230) =] Gaey 9G 
Gy ee 

a Be BOCES os (in) 

q = 176.87 / O,O}0U3/ = 174472 cae im) 
ho = qyy 9, -Fa7 402 2G Wi / rie) 


The following are the results obtained from the computer 


by running the data reduction program (See Appendix F). 


Cie ele oO We ies) 
Ome 7 (6) 
h = 742 (W/m?.K) 


ee 


AP Peo 
UNCERTAINTY ANALYSIS 


Mine same data set (run number 5) that was used for the 
sample calculation was chosen for the uncertainty analysis; 
therefore, the dimensions of the test section, and the meas- 
ured and calculated parameters found in the sample calcula- 
tion were used in this analysis. All uncertainties are 


presented as a percentage of the calculated parameter. 


PeeeeONCERTAINTY IN SOURCE HEAT-TRANSFER RATE 


VI (W) 
I = ie 2,2/ Camp) re 5 ( ame) 


Ve 25V 5 eee ome Oo lias) Ve = es) Oe ae 70) fs) 


A) 
LG 
U1 


uncertainty in measurement and calculation 


On 
I] 


6Q.,/Q, = ((8V,/V,)? + (8I/1)?)?” 


ae Oy - ((0.05/3.58)2 + (0.025/2.27)7)+ 2 
= 1.80 percent 


B. UNCERTAINTY IN SURFACE AREA 


cos) | DOL 

Do = or oon mm.) 6D, = O.1 (mm) 
mee = 203.20 (mm) ol =O eles (min) 
SAL /A, = ((SD,/D,)? + (SL/L)?)*"? 

ony meee (0. 1/15 288 )* > (0.1/203.2)*)*"? 


= 0.63 percent 


eae! 


C. UNCERTAINTY IN WALL SUPERHEAT 


Ate Tegan 


Side 
Ds LO. OO Gaes SE. se iG) 
ho> So 2 22 


tog (D3) / 8 
J Tae 

1 ene nmoceup We reads tee (See Appendix G) 

lav¢ Oe 12 Ae) 


Si eee Yr - a, )4)/7 a ee ee 
mee 


where: 
S.D. = standard deviation 
Since, logarithmic term in eCauasroenmen Ty is too small 
when compared to standard deviation, this term can be 
neglected for uncertainty analysis (i.e., ei TC ee 
T ‘oo ee. 72 Coe ST = 2x $.D. = 5.370 
wo w0O 
= mn di - ° 
ho ee) balers 33.72°= 10 . 00h =e eG) 
6A TAT ((6 BAT)? + Goce 
SAT/AL = ((5.37/23.72)7 > (20s) 7 eee eee 
SAT/AT = 22.7 percent 
D. UNCERTAINTY IN HEAT FLUX 
q = (Q,- 2Q,)/A4, 
oF = 203.17 (WwW) 6Q = 3.65 (W) 
H 
Assuming the same proportion in the uncertanity for Q: 
Q. = 13.15 (W) SO) SOW ) 
F 
Q, - 2Q, = 176.87 (W) 


ar 


s4/q = ( (6Q,AQ- 20.) + (25Q./Q,- 20, )7 + 
(SA /A, )? jee 
Si/q = ( (3665/176-87)* 4 (0.46/176.87)” + 


POO 63°)* «jo 2 
Ze) percent 


Pe iGERTANT TY IN BOLLING HEAT-TRANSFER COEFFICIENT 


Hoe 4) / Ame 


sh / h OG ly ee Gag ly el) ae 


mae) | hh 


ac) 0 eee P0727 )7 
= 22.8 percent 


ia 


Gs 


Tales 
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